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ABSTRACT 


An experimental study was conducted to determine 
the freezing behavior of Ottawa sand and a locally 
occurring silt. The soils were frozen under varied 
pressure, stress history, drainage and temperature boundary 
conditions. 

The results demonstrated that freezing behavior of 
a soil can be predicted by short term freezing tests in 
which heave and net change in sample porewater volume or 
porewater pressure were measured. Results indicate that 
sands expel porewater upon freezing. Freezing behavior 
of fine grained soils is influenced by both soil properties 
and external factors. Freezing behavior of the fine grained 
soils varied from sucking in water leading to ice segre- 
gation and lensing, to expulsion of porewater. It was 
found that the freezing behavior of the fine grained soils 
depended on overburden pressure, stress history, step 
temperature, grain size distribution, porosity and 
permeability. 

Results indicate that heave and heave pressure are 
not a unigue freezing soil property. Experimental results 
support a theory of freezing soils based on the net heat 


flow from a soil sample and on the stress difference across 
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across the ice-water interface. 
Guidelines for a revised frost susceptibility 


criterion are given. 
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CHAPTER ;1L 


INTRODUCTION 


Whenever the mean air temperature drops below 
freezing.for a significant period of time, the effects 
of frost action are of concern in geotechnical engineering. 
To date a large number of facilities constructed in the 
Canadian Arctic have been successfully designed on the 
basis of engineering judgement. Winter construction 
techniques in the more southerly regions have similarly 
been largely developed through experience. With the 
rapidly increasing development of the Arctic regions larger 
and more complex facilities are being proposed for which 
- design precedents are not available and accumulated 
experience is inadequate. There is an urgent need for a 
basic understanding of the behavior of freezing soils. 

Freezing behavior of soils such as heave, heave 
pressure and frost susceptibility, are well known and 
easily identified. The bulk of the work that has been 
carried out on freezing soils is concerned with the 
analysis of the fundamentals of freezing soils. Few 
studies to date are concerned with the phenomenological 
behavior of soils subjected to controlled freezing conditions. 


The purpose of this thesis is to present some observations on 


Sime te ee 
_ : : te ae ” a 


“WOLTOVGOATUI 


woled sqorb sabdaiequsts xis nbom only sowenode 
ajootie add \9mid to boirsq gasottiapse # 302 paisoot? 
-painxnesdipns Levindostosp ai mis9mo5 20.918 nobtos s2eo12 Yo 
sais ot bedauasenos 2oittiios? to xedmun opel 5s otsb of 
ons no fbectetesd viluiessooue nasod ovr pisous m6 Lbeme9 
nofsouztano2 JugntW -tnemsphu pitxesniipas to atesd 
vise t aette syed anotpst yl wendivoe -otom ods at gacdpindoed 
edt daiwW .sonsiisdxs ripyowws bsqolaveb yloprsl-nesd 
uspisl anolpsx sisatA edz to tnemgoleveb patessxoni yibigex 
Hoidw 10% bezoqexq poised sas dole itiost direees: siom bas 
betslumuoos bis sidsfieve ton ste etnebsostq aptesb - 
f tot been taepiy as ei oxoriT .edsupebsat ei sonetueqke 
.eftoa pnissett to xolvsded os to pribastexebau ofasd 
evaor avait 26 dove alioe to roivaded prissext | 
bas awondt [few o16 ,ysilLiditqscveve jeox2 bas sivegstg 
need esi teds Axow oft to ALud SsAT -baitisasbt yLiese 
ads dtiw bamzeonoo 2: alice pnixesti.no tuo beiweep 
wou -eliog paissstt to elejnomsbaut sit io eteylens 
Leviposonsmonedg od? d3 iv beqzeonoo 018 stsb oF esbbude 
_. yattetttbaos patsesxt belloxtio> o3 betoetdve elioa io szoiveried 
7 _ fr anoitsvisedo emos tnseotq oF ei eisedt aid+ to secqiuq odT 


the behavior of freezing soils. 

Initially, the basic outline and goals of the study 
were conceived by McRoberts (1972) as a basis for analyzing 
landslides in Northern Canada. That study was concerned 
with the concept of expulsion of porewater by granular 
soils proposed by Mackay (1972). Field studies of sub- 
surface ice conducted by Mackay (1972) revealed that sand 
and gravel underlay 95% of the ice masses investigated. 
Furthermore, in the remaining 5 per cent of the cases the ice 
was underlain by a clay layer which in turn rested on sands 
and gravels. 

A field case given by Mackay (1972) illustrates the 
porewater expulsion phenomenon. In this instance a large 
lake, which remains partly unfrozen during the winter, is 
partially drained by natural processes resulting in the 
formation of 3 smaller lakes in the original lake basin. In 
the subsequent period of 36 years, 3 pingos developed where 
the smaller lakes occurred. The underlying basin of the 
large lake was also composed of sands and gravels. Mackay 
(1972) believed that as the large lake drained the exposed 
ground surface began to freeze everywhere except under the 
three small lakes. He proposed that the aggrading perma- 
frost resulted in the expulsion of porewater from the lake 
basin sands and gravels resulting in an excess of water 


in areas underlying the as yet unfrozen smaller lakes. 
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However, in time these smaller lakes became sufficiently 
shallow, to freeze completely and to allow the permafrost 
to aggrade in these areas as well. The ensuing result 
was a mound with a core of ice, a pingo. 

The significance of the porewater expulsion concept 
has many important practical implications. Frost heave 
and heave pressure are known to result when water is sucked 
into the soil sample resulting in segregated ice and ice 
lensing. The expulsion of porewater is the reverse effect 
and would indicate that no heave or heave pressure would 
occur during these conditions. The present classification 
of frost susceptibile soils, by Casagrande (1931), is based 
on the relative tendency of a soil to heave or to generate 
heave pressures. If a soil expels water then it would appear 
that it would not heave and hence be classified as non- 
frost susceptible. Clearly the mechanism which controls 
whether a soil will suck in or expel porewater is of 
vital importance in understanding and predicting the behavior 
of freezing soils. 

The following study investigates the concept of 
porewater expulsion by freezing soils. The study concentrates 
on determining which factors, viz., applied pressure, grain 
size, grain size distribution, temperature, and stress 
history, influence the expulsion/sucking in mechanism of 


freezing soils and the relative importance of these factors. 
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Fundamental research into the actual freezing mechanism 
is not part of this study but a mathematical model based 
upon the thermal balance occurring during freezing is 


investigated. 
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CHAPTER) 2 
LITERATURE SURVEY AND FREEZING THEORY 


2.1 Selected Literature Review 


FN tg NES Freezing Mechanism 


The basic mechanism of the freezing of soils is 
well established qualitatively and generally accepted. 
Whenever a negative temperature gradient is imposed on a 
soil, the volumetric heat is removed causing a drop in the 
bulk temperature. When all of the volumetric heat has been 
removed and the bulk temperature is 0 degrees Centigrade 
nucleation of ice occurs with a subsequent release of 
latent heat. If supercooling of the porewater occurs 
nucleation of the ice will occur at some temperature lower 
than 0 degrees Centigrade. The ice crystals grow in the 
direction of heat removal, (Taber, 1929) and impinge and 
apply pressure against any restraining boundary. The pressure 
is relieved by heaving of the soil in the direction of least 
resistance (Kaplar, 1970). The heave is in response to the 
9 per cent volume change at crystallization of the porewater. 
After nucleation the freezing front moves in 
acoordance with the relative balance or imbalance of the 


heat supplied to the heat removed from the soil samole 
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(Anderson and Morgenstern, 1973). If the rate of heat 
removal is greater or less than the rate of heat supplied, 
the freezing front advances or recedes respectively. A 
receding freezing front results in the melting of the 
frozen soil from the freezing front upwards. An advancing 
freezing front results in the engulfment of soil particles 
by ice crystals or in situ freezing. If the rate of heat 
removal equals the rate of heat supplied the freezing front 
becomes immobile, a stationary freezing front. 

Once the freezing front becomes stationary 
and water in sufficient quantity is available to the freezing 
front the ice crystals grow by freezing the adjacent water. 
Due to the rigidity and perfection of the ice crystals, soil 
particles are rejected by the ice resulting in the segregation 
of ice and soil particles (Anderson, 1968). Any change 
however, ne the heat balance/imbalance will result ina 
corresponding change in the behavior of the freezing front. 

Sources of heat to a soil sample are heat conducted 
by the soil sample from the surrounding environment and 
heat released by the porewater, and any water flowing into 
the sample. Heat released by the porewater is volumetric 
heat and latent heat of fusion. The latent heat of 


fusion is*®the*majyor source of “Heat (Nixon, 1972) 
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Basic factors that influence the freezing behavior 
of soils according to Taber (1929) and Penner (1972) are: 
thermal conditions, water supply and nature of the porous 
media. Thermal conditions include the thermal gradient 
imposed on the soil, and the direction or directions of heat 
removal. The amount of heat present and potentially present 
in the soil depends on the degree of saturation, initial 
water content prior to freezing, unfrozen water content 
during freezing, and availability of water in quantity. 
Osler, (1967) expresses the availibility of water under field 
conditions in terms of the depth to the water table from 
the freezing front... If the water table is at a great depth 
and the freezing front is unable to induce water migration, 
closed system freezing conditions are approached. Water 
migration is a necessary condition of ice lensing and 
theories used to explain the mechanisms involved are 
discussed in Section 2.1.B. The nature of the porous media 
depends on the soil type, size and shape of voids, soil 
stress history, and permeability as determined by the grain 
size and grain size distribution. Water solutes and minerals 
present in the soil system also influence the freezing 
behavior of soil. 

The amount of water whose freezing temperature is 
less than 0 degrees Centigrade expressed as a percentage of 


the weight of the soil solids is called the unfrozen water 
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content. Unfrozen water content is an indicator of the 

amount of supercooling occurring as a soil is in the 

process of freezing. Hence depending on the amount of 
supercooling (unfrozen water content) of a soil type the 
freezing may follow or lag behind the 0°C isotherm. 

Neresova (1963) has shown that clean sands freeze at 0 degrees 
Centigrade while clays and silts may freeze at temperatures 
lower than 0 degrees Centigrade. The unfrozen water content 
of a soil is a function of temperature, Anderson (1968), and 

a function of soil type and corresponding surface area 
Anderson (1968), Anderson, Tice and McKim (1973), Dillon and 
Andersland (1966), and Williams (1967). The application of 
pressure at a constant temperature increases the unfrozen 
water content (Hoekstra and Keune, 1967). To date unfrozen 
water contents have been correlated to surface area, 

Anderson and Tice and McKim (1973), Atterberg Limits, freezing 
point depression, clay mineralology, and activity ratio, 
Dillon and Andersland (1966). 

Two equations, the Clausis-Clapeyron and Kelvin 
equations, may be used to describe conditions at the freezing 
front.  Glasstone and Lewis (1964), derive the Clausis- 
Clapeyron equation from basic thermodynamics. The equation 
governs conditions under which a system remains in equilibrium 
at all times as it is being transformed from one phase to 


another as a result of external influences. Changes of 
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state may be caused by the removal of heat or by variations 


in pressure on the system. Anderson and Morgenstern (1973) 


have related the equation to freezing soil conditions in 


the form: 


dP 


oh 


where: 


fs 


change in pressure across the freezing front 


change in temperature across the freezing 


front 


difference in latent heat of the two phases, 


latent heat of fusion of ice 


freezing temperature of the soil 


change in specific volume occurring during 


phase transformation 


The Kelvin equation expresses the pressure relation 


at the interface for small crystals in their own melt. 


Williams (1967), applies the equation to the situation of 


developing ice in a capillary in the form: 
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where: 
PB, = pressure on ice, overburden pressure 
PY = pressure in porewater 
O,, = pressure differential across the ice-water 
interface 
Hears radius of the ice-water mee reoer taken 


positive on the ice side of the interface 


For conditions of optimum heat removal, temperature 
conditions, soil structure and access to a water supply, 
pressures generated by enlarging ice crystals may be 
predicted by the Clausis-Clapeyron equation. Equation 2.1 
may also be used to predict the pressure melting of ice in 
frozen soil from which an estimate of the change in unfrozen 
water content can be made. The Kelvin equation (2.2) can 
be used to predict the decrease (or increase) in the 
porewater pressure whenever the freezing front encounters 
a pore with a radius r, larger than (or smaller than) r._, 


LW 


for a constant overburden pressure. 


Zatix Bane lee Segregation and ice: Lensing 


Segregated ice may be defined as a continuous body 
of pure ice that has developed in but separate from a mass 
of soil particles. Segregated ice may occur in the order 


of microinches in thickness, inches in thickness as ice 
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lenses or up to tens of feet in thickness as massive ice 
bodies. A genetic classification of ground ice has been 
published by Mackay (1972). Segregated ice is generally 
observed in fine-grained soils whether cohesive or non- 
cohesive Taber (1929), Beskow (1935). It has been shown 
that the intensity of segregated ice varies directly with 
initial degree of saturation (Haley, 1953), and varies 
inversely with overburden pressure (Aitken, 1963). 

Ice segregation may occur whenever a continuous 
water supply is available to the freezing front for a 
sufficient length of time under the prevailing temperature 
gradient and pressure differential across the ice-water 
interface. The capability of a soil to supply water to the 
freezing front is a necessary condition for ice lensing. If 
water flow to the freezing front is inadequate, particle 
trapping and ice proliferation may result (Wissa and Martin 
1968). Consequently a lot of work has been carried out 
investigating the mechanism by which a soil is able to 
supply water to the freezing front. Two theories, film 
theory and capillary theory, have been postulated in an 
attempt to explain the phenomenon. The two theories are 
in fact complimentary with the capillary theory being a 
component part of the more general film theory. 

It is now well established that soil particles are 


engulfed by an adsorbed layer of water with different 
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properties than that of the void water, (Anderson. and 
Morgenstern, , 1973) 4. Taber, (1929) was, the first to propose 
this concept while. (Corte,.1962),was the, first,..te,demonstrate 
it-eeathe. thickness of, thesfilmissa function.of temperature 
(Anderson, 1968). 

According to the film theory when a temperature 
gradient is imposed on a soil the water molecules of the 
film crystallize and leave the film. In order to maintain 
the film thickness water molecules from the unfrozen void 
water are absorbed by the film. A process is thereby set 
up to move water molecules to the freezing front. 

Kaplar (1970) postulated the develonment of negative 
film pore pressures. He believed that a quasi-liquid layer, 
‘active film layer’, must be present to provide the drive 
force, for. ce; segregation....The, water molecules: in contact 
with the ice have a high degree of polar orientation with 
the ice lattice and hence have a rigid ordered structure. 
The adsorbed soil water is also formed by strong adsorption 
forces and has a high, structural rigidity. During freezing 
the water molecules in the active film undergo a physical 
readjustment of their position, creating a negative pressure 
ify the, process ,,an order to, fit, into,the,ice orystal, lattice. 
At the time of this readjustment the film force balance is 
disturbed resulting in a rearrangement of the remaining 


film water molecules, and the absorption of void water 
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13 
molecules to restore equilibrium. The demand for replacement 
molecules in the active film layer is governed by the volu- 
metric rate of phase transformation within the layer. 
Freezing of water molecules in the film occurs simultaneously 
with the freezing of void water molecules. 

The initial mechanism in ice lensing is the rejection 
and exclusion of soil particles by the growing ice crystals 
(Anderson and Morgenstern, 1973). Water is drawn to the 
freezing front by negative pressures at the ice-water 
interface and by the imposed negative temperature gradient. 
The ice crystals begin to grow by the freezing of this water 
and=the accumulation of these ice crystals. At crystalli- 
zation the energy that is released is transformed into 
work lifting the overburden, and into heat released at the 
interface which tends to moderate the rate of ice formation. 
The energy involved in the work of the frost heave is due 
to the difference in energy states of ice and water; the 
latent heat of fusion (Anderson, 1968). 

Chronologically the capillary theory was the first 
theory to describe the migration of water to the freezing 
front. The hydraulic theory of Taber (1929) and Beskow 
(1935) was in effect the capillary theory. The basis of 
the theory is that water in small capillaries freezes at 
temperatures lower than 0 degrees Centigrade and that the 


stress difference across the ice-water interface is inversely 
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proportional to the radius of the ice-water interface. 
These relations are given by the Kelvin equation, equation 
2.2, and equations 2.5.01 ven, below. 

In Williams' (1967) development of the theory the 
soil is considered to be a series of capillaries. The 
radius of the capillary (cr) in part controls the behavior 
Ooffethe fneezengerrontinm According sto’ the®capiivary’ theory 
the freezing front continues to advance as long as the 
Capillary radius, ryris greater than the radius of the 
interface, Cay OF X, > Caw When the capillary radius ry 
is smaller than rir the interface radius (rc. < rsh the 
freezing front cannot’enter the capillary’ and its advance 
would bevJhalted. At this®point the cooling ofthe pore 
fluid would continue lowering the fluid temperature below 0 
degrees Centigrade without crystallization. The pore fluid 


would then be supercooled. The depression of the freezing 


point is given by the’ relation: 
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where: 

ea — is the freezing point in degrees Kelvin 
To = normal freezing point 
Nee = specific volume of water at the freezing 


temperature 
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os. = stress difference across the ice-water 
interface 
viw = radius of the ice-water interface taken 


positive on the side of the ice 


L = latent heat of fusion of water 


Due to the reduction in the radius of the ice-water interface 
the porewater pressure would subsequently decrease according 
‘to the Kelvin equatvom, "equatloneyec. (lhis reduction an 
the porewater pressure at the ice-water interface is the 
force that sucks the void water to the freezing point. The 
freezing Of this water results in the growth of ice lénses. 
Jackson and Chalmers (1958, 1966) postulated that 
supercooling supplied the force required to transport water 
to the freezing front. Based on studies of porous media they 
concluded that supercooling controls the free energy of a 
freezing system. Upon the freezing of a soil a reduction 
in free energy occurs in proportion to the amount of 
supercooling. Hence the force that draws water to the 
freezing front is derived from the system free energy 
Change during freezing. 
The capillary theory has several shortcomings and 
is not capable of describing all known freezing conditions. 
Some of these are: (1) It is not universally accepted that 


supercooling is the source of energy that sucks water to 
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the freezing front. Anderson (1968) contends that the 
difference in energy states, the latent heat of fusion, is 

the source of this energy. (2). The capillary theory is 

based on a uniform pore size or uniform particle size 
distribution. Considering soils with a non-uniform particle 
Size distribution, asin virtually all. cases, .difficulties 
arise in assigning a typical capillary pore radius to use 

im the theory. This seriously limits. thesapplicability of 
the theory. (3). The capillary.theory isa static, theory 
based on grain size, pressure and interface stress 
differences... The theory is based on the premise that spore 
radius or grain size and grain size distribution is the 
unigue parameter that controls the ice lensing capability 
of a soil. The theory does not account for any other 
factors affecting the freezing soil system, i.e. temperature 
gradient, and is incapable of accounting for any changes in 


these factors. 


Jule G poke TOs baled Va tg 

When individual ice crystals form and enlarge they 
impinge on confining boundaries generating crystallization 
pressures. Depending on the boundary this pressure is 
often relieved by displacement of the soil surface in the 
direction of least resistance. Displacement usually occurs 


by lifting the overburden or to a lesser extent ina 
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compressible soil by compressing the soil structure. 

As the freezing front penetrates, engulfing the soil 
particles and freezing the porewater in situ, freezing 
pressures and/or heave result entirely from ohase 
transformation. Whenever the freezing front becomes 
stationary however, resulting in ice segregation and lensing, 
heave is a direct result of the buildup of ice lenses. 
Substantial heave as observed in the field can only be 
accounted for in terms of ice lensing and segregation and 
not" as a result of Phases transformation (Taber, 1929).* Ice 
lensing and hence heave may continue indefinitely as long 
as temperature, pressure and water supply conditions remain 
suitable. Whenever one of these conditions is altered, say 
asdecredse inthe water supply, or an increase in temperature 
gradient, the freezing front will advance until a favorable 
balance of conditions is achieved again and the freezing 
front becomes stationary once more. Heaving would be 
minimal while the freezing front is advancing but would 
continue as before once the freezing front stabilizes. 
Repetition of this process results in rhythmic ice banding 
(Maran, «L959 5.2 

The amount and the rate of frost heave is a function 
of many factors. Beskow (1935) and Kaplar (1970) showed that 
heave rate is dependent on the rate of heat extraction or 


freezing front penetration. Penner (1972) found that by 


eaters. Bite | Botany 2 ‘wae 
,gornal a 2a, out LZ) 
éd Ate ees Blot) of ray 


- 
Rh 
7 


7 . 
Saw, rok Apbaues > be © pibens ts peed 


VTeacse>) | Rae KD ‘een Ye - 


Snot om Viesgiditcad suaisacs 1s eon eo 
bees Bio PL Lieiic rhhawe: 157 5w eh Sadat rinte 
yee (Doveagly, gi «Poigchmoo #250) 20° 5ag" ay 

S10 F 64 oqyneet = SeReti. asi ho, iui Srar 
sideaeves € Lisav BoRavbs itiw Jnaez. pai 

SrikesttSd3 Ge Wises Beve.oe Re ona 

Sh bitow Oilvssli paxon sone yishoitasa | 
blugw Jha Britis es 162) ooliesstt ant 
| . San12 primes? ‘apt ginh sxease 
) ’ *; 7 


<2 


18 


increasing the rate of heat extraction, heave rate increases 
to a maximum after which it decreases to a point where 

heave is a result of the expansion of porewater at phase 
transformation. These studies also showed that the response 
of the heave rate to increasing heat removal rates differed 
depending on the soil type. Heaving occurs in the direction 
of heat extraction or crystal growth and total heave was 
reduced by increasing the overburden load (Kaplar, 1970). 
Field studies by Aitken (1963) provided field evidence of 
the effectiveness of surcharge loading in reducing heave. 
Heave .is also dependent on water content (Taber, 1929), 
and.can, -beepravented thypinducing ghigh tensions ei nethe pere= 


water (Gold, 1957),(Penner, 1958, 1966). 


Lag FoR BIS Heaving Pressures 


Upon freezing under constant volume conditions 
pressures are created due to the volumetric expansion 
associated with phase transformation, and due to ice 
segregation and lensing. It is common practice in the 
literature to refer to heaving pressures as those pressures 
which are caused by ice lensing. Heave pressures occur in 
the direction of heat removal and originate at the freezing 
front. Heave pressures are dependent on particle size 
(Gold, 1957), with the smaller particle sizes responsible 
for the maximum heave pressures (Penner, 1968). Miller 


et al. (1960) and Everett and Haynes (1965) developed 
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empirical formulations, based on a study using glass beads 
which related heave pressures to particle or pore size. 

A summary of the findings of several people relating 
Grain size and heave pressure is shown on Figure 2.1. It 


was generally found that the effective diameter, gave 


PAO” 
the best correlation with heaving pressure and so is used 
on the figure. 

The heaving pressures plotted were determined from 
freezing tests with no overburden loading and open drainage 
conditions. Only Sutherland and Gaskin used a back pressure. 
Preparation of the soil samples varied; Hoekstra et al. used 
a compacted sample, Penner densified the sample by vibration, 
Sutherland and Gaskin consolidated the sample from a slurry, 
Yong placed and froze the Sample in a 'loose' condition, and 
Kinosita 'packed' the sample into the soil container. Yong 
also measured heave pressure using proving rings of varied 
stiffness. All samples were reported as being saturated 
prior to freezing. Discrepancies between the data are 
likely a result of less than 100 per cent saturation for 
the duration of the freeze test, variation in sample 
preparation and hence stress history, and a difference in 
freezing rates. Kinosita reported results that were 
determined by both open and closed system freezing tests. 

It would appear however, from the equipment description that 


Kinosita froze under closed drainage conditions. 
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Applying the formulations of Everett and Haynes 
(1965) and the capillary model, Penner (1968) envisaged 
an undulating freezing front which extended as far as 
possible into the larger pores and, which was stopped by the 
smaller pores. He maintained that maximum heave pressure 
was generated when the freezing front everywhere rested on 
the small pores. The tendency for heaving pressures to 
become very large for the smaller particle sizes is shown 
in Figure 2.1. This condition of maximum heave pressure 
likely occurred when the freezing front became stationary 
after a period of rapid penetration. 

Penner (1970) showed that maximum heave pressures 
developed when the soil macrostructure was completely 
broken down and the freezing front was relatively stationary. 
These studies also indicated that the higher the soil density 
the larger the heave pressures. Measurement of heave 
pressures were influenced by the compliance (Kaplar, 1971) 
and small heave movements were sufficient to prevent the 
buildup of large heave pressures in the laboratory. Penner 
(1970) noted that the buildup of large heave pressures in 


the field were dissipated in cracks and discontinuities. 
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22 
2.1.—E. Pore Pressures at the Freezing Front 

The existence of pore pressures at the freezing 
front are clearly demonstrated by the expulsion of pore- 
water (Balduzzi, 1959; Wissa and Martin, 1968), and by 
the attraction of water to the soil sample during open 
system freezing tests. Negative porewater pressures at the 
freezing front are the forces that draw water to the freezing 
front resulting in ice segregation and lensing. 

Based on theoretical developments the magnitude of 
the pore pressures varies inversely with pore size, Gold 
(1957), Penner (1958), and have been evaluated in terms of 
the capillary model by Penner (1966). Williams (1967) 
presented data showing that the air intrusion value could 
be used to estimate pore pressures at the freezing front. 

Williams (1967) froze samples 1 and 2 centimeters in 
height by imposing a temperature of -8 degrees Centigrade 
under open system conditions. Results obtained on 2 natural 
silts with a uniformity coefficient (Cc) of about 80 for 
SrLt. mang C= for- sult 64/15-25, and graded fraction ey ~— ey 
are shown on Table 2.1. 

Sutherland and Gaskin (1973) showed that pore 
pressure at the freezing front in compressible and 
incompressible soils could be reasonably predicted using 
grain size and capillary pore size. Work was done on 


pulverized flyash, coefficient of uniformity C= 3, (referred 
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24 
to as soil 1); Kaolin ot = Suofand: ©7517) isotells #Zieandiz3ssiiand 
a mixture of flyash and kaolin Ce maehis2 pssoisinés. +Hbhe isoil 
specimens were frozen uniaxially by applying cold temperature 
of -17 degrees Centigrade under closed system conditions. 


Results reported show: 


Soil P. P May orOpm in © 
15 Wo 9 W 
KN/M KN/M KN/M 
At 1011 1007 62 
2 1062 1035 465 
3 S015 TOT 491 
4 1080 1080 240 


Data on predicted and measured values of maximum drop in 
porewater pressure for air entry tests were also reported 


as follows: 


Soil Maximum Drop in,Porewater Pressure 
KN/M 
Predicted Measured 
Lt 40 62 
2 100 465 
a 143 491 
4 a, 240 
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Based on these results the authors concluded that predicted 
values of the drop in porewater pressures using the air 
entry method grossly underestimated the actual measured 


Values: 


2olch. Frost Susceptibility ope “Soules 


Fine grained soils and coarse or medium grained 
soils with fines were observed by early workers to be the 
most susceptible to frost action. As a means of classifying 
SOrloein terms some nctimeslsCeDt TDi lstye Omir rOsteactLon ad 
Criterion based on grain size was established. The 
Casagrande criterion for frost susceptibility was one of the 
earliest published. The criterion stated that non-uniform 
soils containing more than 3 per cent of grains smaller 
than 0.02 millimeters and very uniform soils containing 
more than 10 per cent smaller than 0.02 millimeters were 
considered frost susceptible. This criterion was based on 
local materials in New Hampshire and has been widely used 
around the world since 1931. Not all materials in other 
parts of the world, however, behaved similarly to the New 
Hampshire materials and modified criteria were designed to 
match local conditions. As a result a great multitude of 
criteria were developed in various parts of the world backed 
by the experience of people like Beskow, Croney, Schiable, 
Linell and Kaplar, and organizations such as the U. S. Civil 


Aeronautics Administration, and the U.S. Army Corps of 
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Engineers, (Townsend and Csathy, 1963). Each criterion 
met the needs of the person/organization involved but none 
provedijte be ansignif£icant seme eriane over the Casagrande 
criterion and none were as widely accepted. 

All frost susceptibility criteria are based on the 
Capillary model of soil freezing and hence are subject to 
the limitations of the theory. The Casagrande frost 
susceptibility criterion based solely on grain size is 
useful.as a guide to the frost susceptibility of a material 
Hubs ts Inaccurate and incapable of predicting all frost 
action phenomenon for all soil types. An improved frost 


susceptibility criterion is needed. 


4- 2 seelheory /Of Ff reezingsot Soils 


The literature cited above is capable only of 
presenting qualitative and some quantitative theories for 
the formation of ice lenses during the freezing of certain 
soil types. A large gap exists between these theories and 
a comprehensive theory capable of describing the freezing 
of soils as observed in the field and laboratory. A theory 
is presented below in an attempt to fill this gap. 

The theory developed below describes the behavior 
of a soil subjected to freezing conditions. Equilibrium con- 
ditions are assumed to exist when the soil is in quasi-thermal 


equilibrium with the surroundings. This condition occurs in 
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ea) 
a one dimensional system whenever the freezing front is 
stationary with an ice lens actively growing or in the trivial 
Situation of an unfrozen soil with a uniform temperature 
distribution. The theory predicts the behavior of a soil 
when this thermal equilibrium is upset. The theory is 
developed in terms of soil mechanics and heat transfer 


principles. 


2 2iclue Assumptions 


The freezing situation is assumed as shown in 
Figure 2.2. Assumed conditions are: 

(l)} A unit quantity of;soil is being frozen. The 
following soil properties are known or assumed: n, porosity; 


Af water content; Cur coefficient of uniformity Ca DeéQ/Pi9)i 


Cur coefficient of consolidation; k, permeability; Yqr dry 
density, and S, degree of saturation (S = 100%). 

(2) Soil is being uniaxially frozen from the top 
downward with the imposition of a negative temperature at 
the soil surface. 

(3) A net heat balance exists at all times within 
the soil freezing system or Heat out (gq) = Heat in (qi). 

(4) Freezing front elevation is given by the 
distance X, from the top of sample. Soil at a depth less 
than X is frozen while soil at a depth greater than X is 


unfrozen. Unfrozen water content is assumed to equal to 


zero. 
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FREEZING MODEL 


FIGURE 2.2 
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(5) A temperature gradient exists in the soil with 
below freezing temperatures, Ty! at the soil surface, 
freezing temperature, Te, at the freezing front and above 
freezing temperature, Tg? at the base of the soil sample. 
The freezing temperature may be equal to or less than 0 
degrees Centigrade. Heat flow out of the soil sample is 
only due to conduction and is proportional to the temperature 
gradient. 

(6) The water level in the soil is known and 
water flow within the unfrozen soil obeys the Darcy flow 
law: 


3 
q=vA=kia [S31 2.4 


where q is the quantity of water in cubic centimeters 
flowing through the soil; v, is the flow velocity in 
centimeters per sec; A, is the cross-sectional area of the 
soil sample in square centimeters; k, is the permeability 
in centimeters per sec; and i, is the pressure gradient, 
dimensionless. 

(7) Drainage of the soil system is controlled, 
open system or closed system. 

(8) Soil particles are engulfed by a film of 
water that has a state different from the void water. 

(9) A stress difference occurs across the ice- 


water interface. 
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2.2.B. Development of Theory 
Applying surface tension physics to an ice 
lensing system, the pressure difference across the ice- 
water inferface for small crystals in their own melt is 
given by the Kelvin equation, equation 2.2: 


295 
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For the case of freezing soils it appears from the work 
of Williams (1967) and Sutherland and Gaskin (1973) that 
the stress difference P. = ee is nearly constant. For 
the mechanistic model presented here it is postulated 
that the stress difference is equal to or less than a 


certain constant during the growth of an ice lens or 


where K is a characteristic value depending on the soil. 

K may be determined analytically, (Williams 1967; 
Sutherland and Gaskin 1973) or experimentally. To find K 
experimentally two types of freezing tests, open or closed 
system may be run. The closed system tests measures the 
porewater pressure and by substituting this value into 
equation 2.5 K may be determined, On the other hand if one 
assumes that ee in equation 2.5 is equal to zero then K is 
equivalent to the overburden pressure Pie The porewater 


pressure in open system tests is equal to zero whenever no flow 
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of water Pee or out of soil sample occurs. By running 
a series of open system tests to find the stress, Ps, 
at which no water flow occurs, K may be found 

directly. 

Conditions governing the magnitude of a may be 
evaluated in terms of soil mechanics parameters and soil 
system temperature and pressure conditions. 

A basic heat balance may be written for conditions 


shown in Figure 2.2A: 


Heat out = Heat in + Heat of phase change 


or 
ae eg ae p [S34 2.6 
where: 
sia is the heat conducted away from the soil 


sample by the imposition of a negative 


temperature gradient at the soil surface. 


q. - is the heat of the system, volumetric heat 
and any heat conducted into the soil from the 
outside environment by the soil itself or 
by water flowing into the soil. When water 
flows out of the soil a net reduction of qs 


results. 
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In - is the heat of phase change or latent heat 

of fusion that is released by the porewater 

upon oe ee Ora depends on the amount 

of water freezing at the freezing front, 

water content and degree of saturation of the 

unfrozen soil, and the amount of water that is 

sucked into the soil during freezing. q_ can 


p 
be several orders of magnitude larger than qa. 


The heat conducted out of the soil sample is 


proportional to the temperature gradient imposed or 


Io ~ Ae ax A een) =o) 
where K - = CONGCUCTIVity of he frozgen.sail ([cal/°Cem sec] 
dT - temperature gradient that exists in the soil 


sample [°C] 
dX - thickness of soil sample [cm] 
A - cross-sectional area of the soil sample [om] 


The heat of phase change is in direct proportion 


to the volume of water sucked into the soil sample, frozen or 
al - : 
a ¢ (vAt) = Lki A 2.8 
where: 


L - latent heat of fusion of the water flowing to 


the freezing point - [cal/gm] 
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chs) 
V =*volume of®soil frozen [om>] per unit time 


which equals vA = kiA from equation 2.4. 
t - time [sec] 


The pressure gradient, i, is equal to the difference in 
pressure, Ap, causing water to flow to the freezing front, 


through the thickness of the soil sample, z, or: 


From figure 2.2B, Ap equals the difference in pressure 
between the hydrostatic water level at a distance x from the 
freezing front, Pa and the pressure in the water at the 


freezing front, Ee or 
dp =P, - P [gm/cm*] 2.10 


Ee _ 2 
where Py = PL Pp. [gm/em™ ] 2d 


Hence the pressure gradient, i, may be expressed as: 

ee ee BD 

i= — — ele 
Z 
Vw 


Substituting, Equation 2.12 into equation, 2.6: 
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Volumetric heat is proportional to the initial 
soil bulk temperature, the difference between the bulk 


temperature and the ambient air temperature, and the net 
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Specific heat of soil and water flowing into the 


soil is given by: 


Ce ravoLumetricispecific heat of soil 


peal/2c cm? ] given by: 
(0.2 + 1.0 W) 


eis unit weight of dry soil 
3 
[gm/cm™ ] 


W = water content [3%] 


volumetric specific heat of water 


(Gree Cal/ 2G em) 


volume of soil sample ine) 


volume of water sucked into or expelled 


from soil sample [om] 


- initial bulk soil temperature [°C] 


initial temperature of AV fom? ] of water 


neo] 


freezing temperature of soil [°C] 
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Te - freezing temperature of water [°C] 
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Amount of heat, q., conducted by the unfrozen soil sample 


is given by 


pee = A eave 
where 
Ke ~ CONQUCEIVity Of unikrozen soil [cak7°C (cm sec] 


dT - temperature differential existing in the 


unfrozen soil [°C] due to step temperature 


dz - thickness of unfrozen soil [cm] 


A = 


cross-sectional area of soil sample [om] 


Total volumetric heat, Gen is the sum of equations 2.14 
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Substituting equations 2.7, 


equation 2.6 yields: 
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Solving for porewater pressure, Py 
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SuDStI tution of ae determined in equation 2.18 


into equation 2.5: 


enables one to determine whether active ice lensing will 
occur. For example for a Pg =a, ous Ps is greater than K 
(P. > K) no ice lens growth will occur, but if P; is less 
than or equal to K (P; < K) ice lenses could grow. kK is 
determined from subsequent test data. 

The growth of an ice lens may be halted by any one 
of a combination of factors. These factors are: (1) low 
soil permeability, K, resulting in a small insufficient 
discharge to the freezing front; (2) overburden pressure 
high so that the pressure differential across the freezing 
front does not create sufficient suction to induce water 
movement to the freezing front, the pressure differential 
may even be positive resulting in flow away from the freezing 
front; (3) suction at the freezing front may be limited by 
cavitation at minimal overburden pressures such as near 
the surface; and (4) Water supply is inadequate as often 


occurs when the water table is at a great depth. 
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The magnitude of heave, S, may be calculated from 
considerations of a balance of soil porewater volumes. 
Heave is proportional to the increase in sample porewater 
volume upon freezing. Initial sample porewater volume, V., 


is equivalent to (assuming S = 100%): 
V = nv PARA Oe | 


where n isyporosity and V is the total volume of the soil 
sample. Volume of sample porewater after freezing, Vie is 


given by: 


where: 
Ve -~ volume of ice in frozen zone (reduced 


accordingly for unfrozen water) 
V.. - volume of water in unfrozen zone 
V_- volume of ice frozen as a result of water, 


ve being sucked into the soil sample to 


build an ice lens. 


From Figure 2.2 the volume of ice in the soil sample frozen 


to a depth x (assuming all water is frozen) is: 


nv e939 


rh 
>| x 


where 2 is the total depth of sample. Volume of porewater 


in the unfrozen zone is 


nv 2.20 
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38 
where z is the thickness of the unfrozen zone and equal to 


Z= ££ -x. Equation 2.20 then becomes: 


= oa 
Nea ara q any 25 ee 


vee volume of ice frozen as a result of water being sucked 
into tire *so1 P¥to’ buiid™” an’ ice’ Tens“is: 


Vea O09 Fy 2a 
a aw 


where os Bs is volume of water sucked into the soil sample. 
The change in sample porewater volume,AV., may then be 


written: 


AM. (= 4V +t. Vi 


= 109 Fnv+ (1 - %) nv + 1.09 v - nv 
Q aw 
=) O40 SueaatiVies 18 OO; 2s 
Q aw 
AV. may also be written as: 
AVS = SA Pa! 


x 
(0.09 q nv + 1.09 Vine) 225 


Therdepth, of acne freezing 1 ront, x7, aims rogure 2.2 
is proportional to the square root of time. See equation 


4.6. 
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The rate at which a volume of water, Weed is sucked 

into the soil sample is dependent on the rate of heat removal 
and the surcharge loading and stress history of the soil 
at the time of freezing. The volume of water hone may be 
calculated in terms of heat removal. The total heat 
released, Qs, by a volume of water that has been sucked 
into the soil sample is (assuming specific heat of water is 
very small) 
De grata t, = dM, 2.26 


1 


where t - time 
L - latent heat of fusion of water 


From equation 2.6 the heat in, Gis is equivalent to: 


Assuming that qs Of equation 2227 41S the slacentynearror 


fusion then equations 2.26 and 2.27 may be equated to give: 


or 


=o ere 
Vow = o Cae a aes 


The heat out, Qos is given by equation 2.7, and 2 the 
heat of phase change is given by equation 2.13. 
Substituting equations 4.8 and 2.28 into equation 


2.25 heave as a function of time is recovered: 
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The rate of heave, AS, may be deduced from equation 
2.29 by differentiating heave with respect to time: 


ds: anV do 7 ty 
AS = AG ak [0.045 Tr sce esl). re yn 2050 


Total heave and rate of heave are also dependent 
on the freezing temperature of the soil (i.e. unfrozen 
water content). When the soil freezing temperature is 
less than 0 degrees Centigrade the heave and heave rate 


would be reduced accordingly. 
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CHAPTER 3 
TEST PROCEDURE AND DESCRIPTION OF EQUIPMENT 


3.1 Test Procedure 


3.1.A. Materials and Sample Preparation 


Testing was carried out using Ottawa sand and 
Devon silt. A summary of material properties is given 
in Appendix A. Ottawa sand was chosen because of its 
low susceptibility to frost action and in order to verify 
the test program hypothesis. Devon silt was chosen 
because of its high frost susceptibility and general 
similarity to silty materials found in northern regions. 
The bulk of the testing was done on Devon silt. 

Both materials were prepared as a slurry. Devon 
silt was slurried at roughly 1 1/2 times its Liquid Limit 
for ease of handling. The slurries were Geaired to 
ensure saturation, Ottawa sand by vigorously boiling and 
Devon silt by means of a vacuum punp. Distilled water 


was used at all times. 


BASEL Soest. Procedure 
The procedure followed during testing was: (1) 
allow temperature of soil sample to reach a constant 


equilibrium value, (2) apply back pressure to °the soil 
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porewater, (3) impose freezing test boundary conditions, 
andind(4), Ereeze thei soil sample. 

A back pressure was applied to all samples to 
ensure 100 percent saturation. Back pressures of 45 and 
50 pounds per square inch were used. While a back pressure 
was applied, sample temperatures were allowed to stabilize 
and come to equilibrium with the ambient air temperature 
of the cooling chamber or styrofoam cabinet. Samples were 
usually left overnight, for 8,to 14 hours, in.order to 
ensure temperature equilibrium. 

Samples were then consolidated to a desired effective 
stress. All Ottawa sand samples were normally consolidated 
whereas tests were run on both normally consolidated samples 
and overconsolidated samples of Devon silt. 

Prior to freezing the soil sample, desired freeze 
test drainage and displacement conditions were imposed. 

Two drainage boundary conditions were imposed on the soil 
sample, open and closed. Either of two displacement 
boundary conditions, restrained or unrestrained, were further 
imposed on each drainage condition. Open system drainage 
occurred when the soil sample had free access to an outside 
water source. The volume of fluid entering or leaving the 
soil was measured by a volume change indicator, described 
below. During closed system drainage conditions the soil 


sample had no access to an outside water source. Porewater 


cng ler: hnsh i 
ee eer ue : 
oisaeaiy dost athe ea: x ) 
ES omnia eae scutes 
at PeAsquied wis 
ond wstgnae .tensday eon +0 ain gautioon wo 

a sibs a de SF laine te: 


aovidbvette, b4tknab «a ae Solin keane atta sabw wabegea — 
bejebi Lounge vi Liaise Ec asligtee & Wns awedin tie Sebel 
aslani2 beideitomno> vit enter i088 io mud silk ahaa ee , 
dite sovea 40 eslqmse began)! 
essay Bitivsbh \efqmsa Hoe ord trisomnd ot sorr4 ' a 
‘peeogi: sie~ anatdsbads dnemecciqeib Bre ayaniaxh Jaod 
1i6e ats no bewookt oa acho hbneo yrebawed apaniewdk om. 
"select tran 6 Sed 20 sangli .bedels Sas aeqe —— 

tetsu) ataw. beans to benlsutes anit toa 
epankenS mpirays sao r smote eunen speed ‘doe ao: 
shisdye ne a2 ‘cae e622 basi slate Lage 249 nat Sarco | 


ate eatypal 1g Ai asuls to snubov en? 


7 


a 


43 
pressures were measured in this case. Depending upon 
whether the vertical displacement of the freezing piston 
was restrained or not, heave pressure or heave was recorded 
respectively. Methods used to measure pore pressure, 
heave pressure and heave are discussed under the equipment 
section. The bulk of the testing was conducted under 
open system, unrestrained vertical displacement conditions. 

To simulate field conditions, samples were frozen 
uniaxially from the top downward by imposing a negative step 
temperature on the soil surface. The step temperature was 
held constant while the resultant freeze rate depended on 
the temperature gradient and sample water content. A 
constant step temperature was used to facilitate calculation 
and prediction of the advance of the freezing front 
(Appendix B). The movement of the O0°C isotherm was deduced 
from soil temperature data. Other parameters monitored 
were volume change or pore pressure, depending on drainage 


conditions and heave. 


3-2 Description OL Equipment 


3.2.A. General 

The experimental apparatus was basically a modified 
oedometer designed for controlled temperature conditions. 
A schematic diagram of the equipment is shown in Figure 3.1. 


Appendix C lists equipment manufacturers and some specifications. 
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Components of the apparatus were: oedometer, freeze system, 
load and back pressure systems, and the measurement systems. 
The various measurement systems monitored temperature, 
porewater pressure, volume change, vertical displacement, 


and heave pressure. 


3.2.B. Oedometer 

Figure 3.2 shows the oedometer. Noteable features 
are the rolling diaphram, freeze piston, and rubber membrane. 

The rolling diaphram was a rubber impregnated 
figerglas mesh sold commerically under the name Bellofram. 
The diaphram acted as a pressure and water seal around the 
freeze piston, and provided virtually frictionless movement, 
excluding frictional resistance due to the O-ring seal on 
the freeze piston rod, of the freeze piston. The membrane 
used was designed to withstand a 300 pound per square inch 
pressure differential with negligible deformation. The 
placement position of the diaphram provided 1/2 inch upward 
and 2 1/2 inch downward piston travel for a net total of 
3 inches. 

The freeze piston transferred the applied load to 
the soil sample and acted as a heat sink. The sections of 
the freeze piston are shown in Figures 3.3, 3.4, and 3.5. 
The rolling diaphram was cemented to the top of the heat 
exchanger unit. The exchange unit contained the coolant 


circulation maze and functioned as the heat sink. To gain 
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maximum heat transfer efficiency the eleralacian maze was 
located below the rolling diaphram and as close as possible 
to the soil sample; circulation maze length was also maxi- 
mized. 

The freeze piston base plate bore on the soil 
surface and fitted to the bottom of the heat exchange unit 
enclosing the maze. An O-ring was seated around the 
perimeter of the base plate in order to impede the movement 
of the soil slurry upwards past the base plate during 
consolidation loading. Side friction due to the adfreezing 
of this extruded material to the freeze piston was thus 
minimized. To further reduce friction the O-ring was set 
in a slot enabling the O-ring to slide in response to 
movement of the freeze piston. 

The inside of the sample container was lined with a 
thin latex rubber membrane. A thin layer of grease separated 
the container wall and the membrane as shown in Detail 'A', 
Figure 3.2. Due to the inability of the grease layer to 
transfer shear stresses created by the freezing soil, 
friction generated along the sides of the soil sample was 
minimized. Pressures required to extrude the frozen soil 
from the soil container were measured to be as low as 2 


pounds per square inch. 
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S3224C2""Freeze System 


The freeze system is shown in Figure 3.1. Freezing 
of the soil sample was accomplished by imposing a negative 
step, temperature on the soil surface.” This was achieved 
by the circulation of a precooled coolant through the heat 
exchanger maze situated on top of the soil surface. An 
ethylene-glycol and water mixture (50-50 mix by volume) was 
used as coolant. The coolant was circulated by and 
maintained at the desired temperature by a constant 
temperature bath/circulator. The coolant was pumped at a 
rate of 6 3/4 gallons per hour and temperature maintained 
to + 0.1 degrees Centigrade. 

The rate of freezing of the soil sample depended 
on the negative temperature gradient imposed on the soil 


and the rate of circulation of the coolant. 


SE2SDF Loading and Backpressure Systems 


Air pressure, regulated by Nullmatic Pressure 
Regulators and monitored by a pressure transducer, was used 
to apply both load and back pressure to the sample. The 
two air pressures were delivered to the soil sample by two 
separate systems as shown in Figure 3.1. The loading of 
the sample was accomplished by an air pressure chamber 
situated above the freeze piston. Consolidation loads and 
freezing surcharge loads are listed for each test in 


Appendices D, E,F, and G. Application of back pressure to the 
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soil sample was achieved by applying air pressure to a 
water reservoir which was connected to the base of the soil 
container. Back pressures of 40 to 50 pounds per square 
inch, as given on each of the test summaries in Appendices 


D,HEF BP, and Gé 


Re tar 2 Temperature Measurement System 


Thermocouples were used to monitor soil sample 
temperatures. Temperatures were recorded by a Honeywell 
Electronik 15* Strip Chart Recorder. Temperatures were 
measured to + 1/2 degree Fahrenheit. Sample temperature- 
time logs and penetration of O°C isotherm for each test are 
given in Appendices, D, E, F and G. 

Soil sample temperatures were monitored by six 
thermocouples set in the wall of the sample container. See 
Figure 3.2. The thermocouples were vertically spaced l 
centimeter apart with the lowest thermocouple 1.6 centimeters 
from the sample bottom. Starting from the uppermost 
thermocouple the thermocouples were numbered from 1 to 6 
consecutively. Thermocouple temperature data collected 
during each test wereidentified by number in Appendices 
Dye eteeand iGe 

A thermocouple located in the freeze piston base 
plate measured the soil surface temperature. Temperature 
interference due to the base plate was avoided by setting 


the tip of thermocouple below the plat . This thermocouple 
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is labelled To, surface temperature, in Appendices D, E, 
EF} pand ‘G. 

Heat removed from the soil sample was calculated 
from thermocouple data collected in the entry and exit 
lines of the heat exchange unit. One thermocouple, labelled 
Thr was placed in the reservoir of the constant temperature 
bath/circulator and measured the 'cold' temperature of the 


coolant. The second thermocouple, labelled T was situated 


R! 
just past the heat exchanger unit in the return line of the 
coolant circulation tubing. The heat removed from the soil 
sample was in direct proportion to the rise in the coolant 
temperature recorded. 


The cold room and/or styrofoam cabinet temperatures 


were also measured and are shown in Figure 3.6. 


Sore Las Pore Pressure Measurement system 


Porewater pressures were monitored by a pressure 
transducer located in the base plate of the sample container 
as shown in Figure 3.2. The transducer was calibrated to 0.01 
pounds per square inch. Pore pressure measurements were 


taken during closed system freeze tests. 


3.2.G. Volume Change Measurement System 


Volume change of the soil sample, in terms of pore 
fluid entering or leavinc, was measured by a 25 cubic 


centimeter burette. Back pressure could be applied to 
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55 
the burette system, as shown in Figure 3.1. Volume change 


measurements were taken during open system freeze tests. 


Se enue. Vertical Displacement Measurement System 


A dial gage measured the vertical displacement of 
the soil surface to 0.001 inches. Soil surface displace- 
ments were monitored during consolidation and freezing tests 


for both open and closed system tests. 


3.2.1. Heave Pressure Measurement System 


A load cell was designed to attach to the freeze 
piston rod and to be fixed to a rigid frame to measure 
anticipated heave pressures. The load cell was made of 
machined aluminum and calibrated to measure + 0.01 pounds 


per square inch. The load cell is shown in Figure 3.1. 


Shears Equipment Assessment and Recommendations 


The rolling diaphram performed satisfactorily and 
added a high degree of flexibility to the type of freeze 
testetnat, could be conducted.) initial’ familiarization, 
installation and sealing of the diaphram proved quite 
troublesome. The system used to anchor the diaphram to 
the freeze piston could be improved to provide faster 
installation of new diaphrams. Under the present arrange- 
ment with the diaphram cemented to the freeze piston, a 
day is required to replace a damaged diaphram, and time 


is lost. In the original design the bellofram was anchored 
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56 
by fitting the freeze piston base plate and heat exchange 
unit together, squeezing the bellofram between the two 
sections. This system proved unsatisfactory and was 
abandoned. This problem is minor in nature but to reduce 
lost experimental time the latter system is more desirable. 

The attainment of a 100 per cent pore pressure 
response (B value equals 1.0) to the applied consolidation 
loading proved to be a problem and an area that should be 
improved. The crux of the problem concerns the rolling - 
O-ring seal on the freeze piston base plate (Figure 3.5). 
The O-ring was set to provide maximum impedence to the soil 
slurry from moving upwards past the base plate during 
consolidation loading without creating excessive normal 
pressures and hence friction against the sample container. 
The O-ring seating used satisfactorily impeded drainage and 
generated acceptable amounts of friction but was incapable 
of confining the pressure within the soil chamber. 
Consequently the air space between the rolling diaphram and 
freeze piston base place introduced a compressible component 
into the pressurized zone resulting in a finite reduction 
of the soil container pressure. This problem occurred only 
during the consolidation loading stage and ceased to exist 
once a plug of frozen soil develops. A new design is needed 
that will confine the soil pressures without creating 


excessive side friction. One solution might be to attach 
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57 
the rubber latex membrane surrounding the soil sample to 
the freeze piston. 

Rate of heat removal is another area that could 
be improved. The heat removal capacity of the freeze 
system is a direct function of the rate of coolant 
circulation. Higher flow rates could be achieved by 
reducing the resistance to flow by using large diameter 
circulation tubing and heat exchange maze flow channels, 
and by using a less complex circulation maze in the heat 
exchange unit. These diameters must be optimized with 


respect to dimensions of the freezing piston. 
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CHAPTER 4 
TEST RESULTS AND CALCULATIONS 


4.1 Introduction 

A general outline of testing procedure was given 
in Chapter 3. Details of the testing program, and 
processing of the data are given below. 

Data for each soil type are given in the following 
appendices: Ottawa sand data are summarized in Appendix D, 
Devon silt data are listed in Appendices E and F, and 
modified Devon silt data are summarized in Appendix G. 
Whenever more than one type of test was conducted on a soil 
type, the testwas identified by the series number and test 
number. For example, Test E-1-2 refers to test 2 of series 
1 listed in Appendix E. 

Results of calculations are also summarized below 


with typical calculations given in Appendix B. 
4.2 Test Results 


4.2.A. Ottawa Sand 

Three freezing tests were run on Ottawa sand to 
verify the concept that during freezing cohesionless 
granular soils expel water. Soil properties of Ottawa 


sand are given in Appendix A. Summaries of soil data and 
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freezing test data are given in Tables 1 and 2, Appendix 
D, respectively. A summary of freezing test results for 
Ottawa sand is given in Figure 4.1. 

All freezing tests carried out on Ottawa sand 
were conducted under open drainage and unrestrained heave 
conditions. Details of stress history and freezing step 


temperature for each test are as follows: 


Test. D=1 

A slurry was placed in the soil container and 
allowed to consolidate under an effective stress of 1.04 
Kg/om*. A step temperature of -10°C was imposed on the 
soil surface. Freezing test data are given in Appendix D, 


(Figure D-1l). 


TestsD—2 
A slurry was placed in the soil container and 
densified by tapping the sides of the container. The 


purpose of this test was to investigate the influence of 
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porosity on the freezing behavior of sand. The slurry was 


allowed to consolidate under an effective stress of 1.04 
Kg/om*. A step temperature of -10°C was applied to the 
soil surface. Data for this test are given in Appendix 


D, ((Figure-—D-2).. 
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LlestaD-S 
A slurry was placed in the soil container and 
allowed to consolidate under an effective stress of 0.50 
Kg/om. This test was run to check the influence of 
overburden load on freezing behavior. The soil was s 
subjected to a -9.39°C step temperature. Freezing test 


data are given in Appendix D, (Figure D-3). 


At the end of these tests the dimensions of the 
sample container and arrangement of thermocouples in the 
container wall were changed to those shown in Figure 3.2. 
The cross-sectional area of the sample container was 
changed from 81 cm”, used above, to 83.4 cm? . The height 
of the sample container was reduced to accommodate the freeze 
piston and the large consolidation strains of the silts 
which were used for the remainder of the testing program. 
The thermocouples were arranged in three sets of pairs in 
the above tests. The thermocouple arrangement shown in 
Figure 3.2 was adopted to gain a more detailed soil 
temperature profile during freezing. Thermocouple Tee 
was also installed in the return line of the glycol 


circulation tubing to monitor the amount of heat removed. 


4.2.B. Devon Silt 
Devon silt was used to test the porewater expulsion 


concept on fine grained soils. A summary of the mechanical 
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properties: of this soil is given in Appendix A.° Freezing 
tests were conducted on normally consolidated and over- 
consolidated samples. Summaries of data obtained for 
normally consolidated samples are given in Appendix E, 
(Tables 1 and 2), and for overconsolidated samples in 
AppendixvF, Tables (1 and 2). Void ratio versus effective 
stress plots are given in Appendix E, (Figure E-0-1) for 
normally consolidated samples and in Appendix F, (Figures 
F-0-l, F-0-2, F-0-3) for overconsolidated samples. Details 


of test series run on Devon silt are as follows: 


(1) Test Series E-1l 

The purpose of Test Series E-1l was to investigate 
the freezing behavior of normally consolidated Devon silt. 
Freezing tests were run at varied effective stresses and 
void ratios. As summarized in Appendix E, (Table 1) a series 
of freezing test, Tests E-l-1 to E-1-7, were carried out 
under effective stresses ranging from 0.72 to 2.42 Kg/om@. 
Freezing test data for each of these tests are summarized 
im Appendix: Bynd(Tablei:2)\.. WA ‘summary? of ‘test? results for 
this test series is given in Figure 4.1. 

For each test of this series a slurry having 
a liquidity index between 3.5 and 4.0 was placed in the 


sample container and consolidated to the desired effective 


stress. Consolidation data are summarized in Appendix E, 
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(Table 1 and Figure E-O-l). The samples were frozen by 
imposing a step temperature of -4.44°C on the soil surface. 
All freezing tests were conducted under open drainage and 
unrestrained heave conditions. Data for tests E-1l-1l to 
E-1-7 are shown in Appendix E, (Figures E-1-1l to E-1-7 
respectively), 

Test E-1-2 was conducted at roughly the same 
effective stress as Test E-l-1l to demonstrate repeatability 


of results. 


(2) Test Series E2 

This test series was conducted under closed 
drainage conditions in order to investigate the freezing 
behavior of normally consolidated Devon silt in terms of 
porewater pressure generation. It was anticipated that 
the freezing behavior mapped in Test Series E-l and E-2 
would be consistent. 

In a procedure similar to Test Series E-l, 
freezing tests E-2-1 to E-2-5 inclusive of series E-2 were 
carried out under effective stresses ranging from 0.97 to 
1.47 Kg/om respectively. For each test a slurry having a 
liquidity index between 3.1 to 3.9 was placed in the sample 
container and consolidated to the desired effective stress. 
Consolidation data are summarized in Appendix E, (Table l 
and Figure E-O-l). A step temperature of -4.44°C was applied 


to the soil surface and the soil frozen under unrestrained 
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64 
heave and closed drainage conditions. Data for this series, 
Tests E-2-1 to E-2-5, are shown in Appendix E, (Figures 
E-2-1 to E-2-5 respectively). A summary of test results for 
this test series is given in Figure 4.2. 

Tests E-2-4 and E-2-5 were run at roughly the 
same effective stress to demonstrate reproducibility. Pore- 
water pressure results of Test E-2-2 cannot be fully 


explained. 


(3) Test Series E-3 

Boundary conditions of both tests in this series 
were changed while the freezing process was in progress. 
Freezing was initiated in Tests E-3-1 and E-3-2 under open 
drainage and unrestrained heave conditions and were continued 
until a consistent set of porewater expulsion data points 
was generated. At this time the sample drainage was stopped 
and the freezing test continued under closed drainage 
conditions. The purpose of these tests was to attempt to 
Measure the porewater pressure (closed drainage conditions) 
that caused the porewater expulsion under open drainage 
conditions, and to see if freezing conditions were different 
for open and closed drainage cases. 

Again for each test a slurry at a liquidity index 
of 3.8 and 3.5, for Tests E-3-l1 and E-3-2 respectively, was 
placed in the sample container. The slurries were then 
consolidated under the chosen effective stresses. The 


consolidation stress of 1.47 Kg/em* for Test E-3-1 and 1.96 
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Kg/em* for Test E-3-2 were chosen to ensure expulsion of 
porewater during freezing. (These values were based on 
results obtained in Test Series 1). Consolidation data 
are given in Appendix E, (Table 1 and Figure E-O-1). 

As stated above freezing was initiated in both 
tests under open drainage and unrestrained heave conditions, 
and continued until the porewater expulsion rate became 
constant. The drainage was then closed in both tests and 
porewater pressures were measured. Once consistent pressures 
were being recorded the drainage system was opened in both 
tests, allowing free drainage. Test E-3-1 was terminated at 
this stage. After a period of free drainage the drainage 
system was closed for a second time in Test E-3-2 and pore- 
water pressures measured once again. Step temperatures of 
“4.44 Geana —o.00°G were usedoin Tests B=3—)7and &—-3-2 
respectively. Data collected for these tests are given in 


Appendix E, (Figures E-3-1 and E-3-2). 


(4) Test Series F-l 

In Test Series F-l a sample of Devon silt was 
overconsolidated to investigate the influence of stress 
history on freezing behavior. Series F-1l is the first of 
several test series conducted to investigate this relation- 
Snip. 

A slurry at a liquidity index of 3.6 was placed 


in the sample container and consolidated under an effective 
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stress ofe7&70 Kg/em*. The applied stress was then reduced 
Ome. Kg/em* and the sample allowed to rebound. A number 
of freezing tests were subsequently conducted at various 
effective stress conditions. All freezing tests. of this 
series were run under open drainage and unrestrained heave 
boundary conditions. Details of stress path followed and 
applied step temperature for each test of the series are 


as follows: 


Testarb=l=i 

Prior to freezing in this test the sample was 
allowed to consolidate under an effective stress of 7.70 
Kg/em* and then allowed to rebound under an effective 
stress of 0.1 Kg/em? as described above. Consolidation 
data are presented in Appendix F, (Figure F-O-1). A step 
temperature of -4.44°C was applied to the soil surface. 
The data for this test is presented in Appendix F, (Figure 
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Afterscompletionlofetest¥F=1=l1.the effective stress 
was increased to 0.25 Kg/cm? and the sample allowed to 
consolidate under the higher load. Consolidation data are 
summarized’ in Appendix F, (Figure F-0-1). A surface 
temperature of -4.94°C was applied to the soil surface. 
Test data for Test F-1-2 are shown in Appendix F, (Figure 
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TestsP=-l=3 
The applied effective stress was doubled to 0.5 
Kg/cm* and the sample allowed to consolidate. Consolidation 
data are presented in Appendix F, (Figure F-0-l). A step 
temperature of -4.44°C was imposed on the soil surface. 
The data for this test are presented in Appendix F, 


(Figure F-1-3). 


Test F-1-4 

The applied stress was increased to 0.74 Kg/em* and 
the sample allowed to consolidate. A step temperature of 
-5.00°C was imposed on the soil surface. This test was the 
first to show that porewater could be expelled by freezing 


an overconsolidated silt. Data for this test is contained 


in Appendix F, (Figure F-1-4). 


Testes 

The applied effective stress was increased to 0.98 
Kg/em?, the sample was allowed to consolidate and then 
frozen. A step temperature of -4.79°C was used. Test data 


is given in Appendix F, (Figure F-1-5). 


A summary of test results for this test series is 


given in Figure 4.3. 
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(oye Dest sock les! 2s 

In this test series the sample used in Test 
Series F-l was frozen under closed system drainage 
conditions. The purpose of this series was to determine 
whether the porewater expulsion phenomenon of Test Series 
F-1 could be duplicated in terms of porewater pressures. 

Consolidation data are shown in Appendix F, 
(Figure F-0-1). All tests of series F-2 were conducted 
under closed drainage and unrestrained heave conditions. 
Stress path and surface temperature details for each test 


are as follows: 


Tes telw2—. 

Prior to freezing the applied effective stress was 
reduced from 0.98 Kg/em? of Test F-1-5 to 0.50 Kg/em* and 
the sample was allowed to rebound. The sample was frozen 
by imposing a step temperature of -4.44°C on the soil 
surface. Test data are presented in Appendix F, (Figure 
F-2-1)% 

testelh=2=2 

The applied effective stress was increased to 0.98 
Kg/em? in an attempt to generate positive porewater pressures 
during freezing. After consolidation was complete a step 
temperature of -4.94°C was applied to the soil surface. 


Freezing test data are given in Appendix F, (Figure F-2-2). 


Test F-2-3 


The applied effective stress was reduced to 0.70 
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Kg/cm”, the sample was allowed to rebound, in an attempt 
to repeat the generation of negative porewater pressures 
measured in test F-2-2. A step temperature of -5.00°C 
was used. Data for this test are given in Appendix F, 


(Figure F-2-3). 


A summary of test results for this test series is 


given in Figure 4.2. 


(6) Test Series F-3 

For this test series the step temperature used 
in previous test series was roughly doubled. The sample 
used previously in Test Series F-1 and F-2 was also used 
for this series. Consolidation data are given in Appendix 
F, (Figure F-0-1). All freezing tests were conducted under 
open drainage and unrestrained heave conditions. Stress 


path and surface temperature data for each test are as follows: 


Test Pe3eu. 

The effective stress of 0.70 Kg/om used in Test 
F-2-3 was maintained and used in this test. The sample was 
frozen by applying a step temperature of -9.72°C. Test data 
areéegiven in-Appendix)\F,. (Figure .F-3-1). Due to,the reduced 
time required to freeze the soil sample all summarized data 
were extrapolated to time t = 160 minutes for comparison 
with previous data. 


Test F-3-2 


The effective stress was reduced to 0.25 Kg/om 
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and the soil allowed to rebound. A step temperature of 
-9.72°C was applied to the soil surface. The data for this 


test is presented in Appendix F, (Figure F-3-2). 


Test F-3-3 

The soil sample was consolidated under an 
effective stress of 0.50 Kg/em*. The soil surface was 
subjected to a -9.78°C step temperature. Data for this 


test is presented in Appendix F, (Figure F-3-3). 


A summary of test results for this test series 


is given in Figure 4.4. 


(7) Test Series F-4 

In this series freezing tests were run ona 
moderately overconsolidated Devon Silt sample. The 
purpose of the series was to further investigate the 
influence of stress history on freezing behavior. 

Prior to freezing the soil sample was 
consolidated under an applied effective stress of 4.88 
Kg/em? and then allowed to rebound under a stress of 0.46 
Kg/cm*. Consolidation data are summarized in Appendix F, 
(Figure F-0-2). All freezing tests in this series were 
conducted under open drainage and unrestrained heave 
conditions. Details of stress path and surface temperature 


for each test are as follows: 
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Test F-4-1 

A slurry having a liquidity index of 3.6 was placed 
in the sample container, allowed to consolidate then rebound 
under effective stresses of 4.88 and 0.46 Kg/cm* 
-respectively. The sample was subjected to a step temperature 
of -7.36°C for the first 30 minutes and of -4.83°C for the 
remainder of the test. The change in the step temperature 
was a result of problems experienced with the regulation 
of glycol flow prior to the test. Due to the collapse 
of the circulation tubing the flow rate was significantly 
reduced. In order to maintain the same rate of heat 
extraction used in previous tests the step temperature was 
increased in proportion to the decrease in glycol flow. 
However, as test F-4-1 progressed it became evident, from 
frost penetration data, that a glycol flow rate roughly 
equivalent to earlier tests had been achieved by the 
equipment modifications. The step temperature was then 
increased to -4.83°C. Data for this test are given in 


Appendix F, (Figure: F-4=1). 


Test F-4-2 

The applied effective stress on the sample was 
increased to 0.82 Kg/em? and the sample was allowed to 
consolidate. The sample was frozen by applying a step 
temperature of -5.24°C. Data for this’ test is presented 


in Appendix Fi, (Figure F-4=2). 
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Test F-4-3 

The effective stress was increased to 1.07 
Kg /om* in an attempt to cause porewater expulsion during 
freezing. The sample was allowed to consolidate and then 
frozen by a step temperature of -4.72°C. After roughly 
50 minutes of testing a city wide power failure occurred 
and the test was necessarily terminated. The soil sample 
was removed because of feared contamination by kerosene 
from the volume change indicator. Because of the 
consistency of the expulsion porewater results and the 
favorable agreement with earlier tests, the author 
extrapolated the data to time t = 160 minutes. The data 


for this test is presented in Appendix F, (Figure F-4-3). 


Because the results of Series F-4,tests F-4-1 to 
F-4-3, and as a whole were consistent with previous 
observations this test series was not repeated. A summary 


of. test results. for this test.series’/ is given in: Figure 4.38. 


(8) Test Serres F=5 

In Test Series F-5 a Devon silt sample was 
heavily overconsolidated, prior to freezing, to complete 
the study of stress history on freezing behavior. Heave 
was monitored in the first 4 tests of this series; heave 


pressure was measured in the last test, Test F-5-5. 
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A slurry at a liquidity index of 3.7 was placed 
in the soil container and allowed to consolidate under an 
effective stress of 8.79 Kg/com*. The applied stress was 
then “reduced to 0.1 Kg/em* and the sample allowed to 
rebound. Consolidation data are shown in Appendix F, 
(Figure F-0-3). Details of stress path and freezing 


boundary conditions for each test are as follows: 


Test F-5-1 

As mentioned above the slurry was consolidated then 
allowed to rebound under effective stresses of 8.79 and 0.1 
Kg/om* respectively. A step temperature of -5.56°C was 
applied to the soil surface and the sample frozen under 
open drainage, unrestrained heave conditions. Data for 


this test are,given in Appendix F, (Figure F-5-1l). 


Test F-5-2 

The applied effective stress was increased to 0.24 
Kg/cm?, the sample was consolidated and then frozen under 
open drainage and unrestrained heave boundary conditions. 
A step temperature of -4.78°C was used. Test data are 


shownhineAppendixeE, (Figurne.<F-5-2). 


Testeu-o-3 
The sample was allowed to consolidate under an 
effective stress of 0.49 Kg/cem-. The sample was subjected 


to a step temperature of -5.00°C and frozen under open 
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drainage and unrestrained heave conditions. Data for this 


test are given in Appendix F, (Figure F-5-3). 


Test F-5—4 

The applied effective stress was increased to 0.97 
Kg/em* and the sample allowed to consolidate. A step 
temperature of -5.17°C was applied to the soil surface and 
the sample frozen under open drainage and unrestrained 
heave conditions. Test data are shown in Appendix F, 


(Figure F-5-3). 


A summary of test results for tests F-5-1 to F-5-4 


of this series is given in Figure 4.3. 


Test. h—o=5 

In this test heave pressure was measured by 
restraining the heave displacement of the soil surface 
using a load cell fixed to a rigid load frame as shown 
in-Pigure 3.1. 

The purpose of this test was to measure heave 
pressure of a sample at an effective stress at which no 
net change in sample porewater volume occurs. Based on 
results of previous tests of this series an effective stress 
of 0.69 Kg/em? was applied to the sample. The sample was 
allowed to rebound under this reduced stress and frozen under 


open drainage and restrained heave conditions. A step 
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temperature of -5.17°C was used. Any vertical movement of 
the freeze piston due to load frame compliance was monitored 
for control purposes. Test data are shown in Appendix F, 


(Figure F-5-5). 


4.2.C. Modified Devon Silt 

The grain size distribution of Devon silt was 
artifically changed in order to study the influence of grain 
size distribution on freezing behavior. Modification of 
the characteristic grain size distribution was accomplished 
by increasing the proportion of finer silt sizes and clay 
sizes. The finer particle sizes were removed from a second 
Devon silt sample by elutriation. Grain size distribution 
plots of typical Devon silt and Modified Devon silt are 
shown in Appendix A, (Figure A-1). 

A slurry prepared at a liquidity index of 5.1 was 
placed in the soil container and allowed to consolidate 
and rebound under effective stresses of 7.7 and 0.1 Kg /em* 
respectively. Consolidation data are contained in Appendix 
G, (Table 1 and Figure G-0). Details of stress path and 


freezing boundary conditions for each test are as follows: 


Test G-l 
As described above the slurry was allowed to 
consolidate and rebound under effective stresses of 7.70 


and 0.1 Kg/em? respectively. A step temperature of -5.22°C 
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was applied to the sample surface. Freezing occurred under 
Open drainage and unrestrained heave conditions. Data for 


this test are given in Appendix G,.(Figure-G-1l) ; 


Test Gaz 
The effective stress was increased to 0.24 Kg/em* 
and the sample allowed to consolidate. Freezing was carried 
out under open drainage and unrestrained heave conditions. 
A step temperature of -4.87°C was used. Data for this test 
are given in Appendix G, (Figure G-2). 

The volume of porewater sucked in during this test 
was less than anticipated. The validity of the test data 


is open to question due to a suspected leak in the volume 


change indicator. 


Test G-3 

The effective stress was increased to 0.49 Kg /em? 
After sample consolidation was completed the sample was 
frozen under open drainage and unrestrained heave conditions. 
A step temperature of -4.72°C was used. Data for this test 
are given in Appendix G, (Figure G-3). The porewater 
volume results for this test are also questionable as the 
leak in the volume change indicator was not successfully 


stopped. 
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Test G-4 

The sample was consolidated under an effective 
Sshressacte0a73 Kg/cem-. Step temperature of -4.79°C was 
applied to the soil surface. The soil was frozen under 
Open drainage and unrestrained heave conditions. Test 
data are contained in Appendix G, (Figure G-4). The pore - 
water volume results are also questionable for this test. 
Again the leak in the volume change indicator was not 


successfully stopped. 


Test G-5 

The effective stress was increased to 0.98 Kg/em* 
and the sample allowed to consolidate. Freezing was 
conducted under open drainage and unrestrained heave 


conditions. A step temperature of -4.89°C was used. Data 


for this test are shown in Appendix G, (Figure G-5). 


Test G-6 

The effective stress was increased to 1.22 Kg/em? 
and the sample allowed to consolidate. Freezing was 
conducted under open drainage and unrestrained heave 
conditions. A step temperature of -5.00°C was used. Data 
(ope “itil test are given in Appendix G, (Figure G-6). Results 


of this test are consistent with Tests G-l, and G-5. 
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Test G-7 

The effective stress was increased to 2.23 Kg/em* 
in an attempt to expel porewater during freezing. The 
sample was allowed to consolidate under this pressure 
and frozen under open drainage and unrestrained heave 


conditions. A step temperature of -5.00°C was used. Data 


for this test are shown in Appendix G, (Figure G-7). 


A summary of test results for tests G-l to G-7 


of this series is given in Figure 4.1. 


Test G-8 

Test G-8 consisted of measuring the heave 
pressure generated by a freezing soil with free drainage. 
As in Test F-5-5 the effective pressure at which this test 
was carried out at would not result in any net change in 
sample porewater volume. Based on results of Tests G-l, 
G-5, 6, and 7, an effective stress of 2.24 Kg/cm” was used. 

The sample was allowed to rebound under this 
effective stress of 2.24 Kg/em*. The sample was subjected 
to a step temperature of -4.50°C and frozen under open 
drainage and restrained heave conditions. Vertical movement 
of the freeze piston, due to compliance, was monitored for 
Control purposes. Test data for this test are contained in 


Appendix G, (Figure G-8). 
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The apparent correlation of results of Tests G-l, 
2, 3 and 4 is recognized by the author. The leak that 
occurred in the volume change indicator during Tests G-2, 
3 and 4 was very small but its possible effect on test 
results must be acknowledged. Because of this possibility 
of error and the satisfactory duplication of results of 
Tests G-5 and G-6 the relationship shown in Pigure 4.1 is 


recommended. 


4.3 Calculations 
Soil property and freezing test calculations are 


presented below. 


4.3.A. Soil Property Calculations 


Initial height of the slurry sample was the average 
of nine vernier measurements. The height of sample at any 
time during testing was determined by adding or subtracting 
the net change in the dial gage reading to the initial 
height. 


The void ratio, e, at any time is: 


oO 


H 
o 


where H is the sample height [cm] 


H is the height of soil solids [cm] or 
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where W., rspthetweightwof therdrytsoil “['cem] 


Geis thegspeciticlgravity of the soil solids 


line is the unit weight of water = L[gm/cm?] 


A is the cross-sectional area of the sample 


container heme (A’= 81 aa for sand and 


A = 83.4 em* for the silts). 


The degree of saturation at the beginning of the 
test, S. was calculated from the initial water content and 


void ratio by 


WS 
So af ie Ap 
re) 
where ihe is the initial water content 


ey is the initial void ratio. 


oe for virtually all tests was equal to 100%. 


The coefficient of consolidation, Coe was determined 


from e - log time plots in the conventional manner: 
0.848 ie 2. 
Cc = [cm”/sec] a 
Vv t 99 


where H is the average sample height 


too is the time at which 90% of total consolidation 


has been completed. 
Actual t 99 data are of little significance to this study 
and have been omitted from data summaries. 
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Permeability, k, was then determined by: 
= ye Cc mM [cm/sec] 4.4 


where m is the compressibility 
i 2 
WEE prog” Hest [cm™ /gm] 4.5 


Prior to the freezing tests, void ratio ee, was 
Calculated from equation 4.1 by substituting the sample 
height prior to freezing for H. Porosity was then 


calculated: 


As shown in Table 1 in each of Appendices D, E, 
Fland:G, So calculated from equation 4.2 was virtually 
equal to 100 per cent in all cases. Consequently water 
content prior to freezing, Wer was calculated from equation 


4.2, assuming So = 100% or 


x 100% 4.7 


e 
ut G 


Ae Be Hreeging Test Calculations 


(1) Penetration of 0°C Isotherm 

The movement of the freezing front in a soil 
subjected to a step decrease in temperature has been 
described by Neumann and is given by Carslaw and Jaeger 


(1947). Assuming that the properties of the frozen and 
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thawed regions are homogeneous and independent of 


temperature, the movement of the 0°C isotherm through the 


unfrozen soil is: 


where 


where 


xX = a ft 4.8 


X is the depth of frost penetration 
t@ais time 
a is a constant which is determined as a root of 


the transcendental equation 
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k:, k. are the diffusivities of unfrozen and frozen 
soil fom~/sec} 
Ku K.. are the thermal conductivities, of unfrozen 


and £rozen soil | Ccaly °C) cm sec] 


Cay, Ce are the volumetric heat capacities of 
unfrozen and frozen soil [cal/°C ane soil] 


te is the uniform initial ground temperature [°C] 


is is the applied constant surface step temperature 
[°c], L is the volumetric latent heat of the 
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Assuming a linear temperature distribution in the 
frozen zone and ignoring the temperature profile in the 


unfrozen zone Stefan solved equation 4.8 in the form 


okt 
X = ——— fe acolo 


from which the constant term can also be written in the 


dimensionless form: 


st 
O, = 


IK 2 
u 
Le 


Cer 
= ; S and is called the Stefan number. 


2 
where Ss 


Volumetric heat capacity of unfrozen soil, Cul and 


volumetric latent heat of the soil, L, may be determined as 


follows: 
1.0w 
Ch F ota (C + i60 ) 4.12 
where 1.0 = heat capacity of water [cal/gm°C] 
— t rt 
L Tonle Amat Ww) L 4.13 


The terms of equation 4.12 and 4.13 are defined as: 


Cn denotes the heat capacity of the soil grains 
[calL/gm? Cc] 


Yq denotes the dry density of the soil [gm/om?] 


w denotes water content [$%] 
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ie denotes unfrozen water content [gm water/gm 
(ice + water) ] 


L' denotes latent heat of water [= 79.6 cal/gm] 


The test equipment was designed to impose a step 
temperature on the soil surface. Figure 4.5 shows the 
penetration of the 0°C isotherm for several selected freezing 
tests. Penetration data of Figure 4.5 was monitored by 
thermocouples. Appendix B, (Table 1) summarizes the X, vt 
data plotted in Figure 4.5. 

The curves shown in Figure 4.5 are bi - and 
trilinear and hence do not follow the relation described by 
equation 4.8. An average rate of penetration, Qos is 
calculated for each selected test in Appendix B, (Table l). 
Table 4.1 summarizes the calculated a, values. 

For each test shown in Figure 4.5, theoretical a 
values were calculated using equation 4.10. Unfrozen water 
contents we) of 0 and 5% were assumed in the calculations. 
A sample calculation is contained in Appendix B, (Section 
1). A summary of calculated a values is given in Table 4.1. 

To this point, the discussion has been concerned 
with depth of frost penetration as determined by thermo- 
couple measurements of the sample temperature distribution. 
For normally consolidated tests the penetration of the 


freezing front was also measured by Vernier at the end of 
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OS - Ottawa Sand 


TABLE 4.1 


Rates of Penetration of Freezing Front 


Test 


Experimental 
Os [cm/ min] 
Thermocouple Measured 

= 05323 

= 0.341 
0.342 0.267 
0.407 OAH) 
0.376 O27 
0.384 = 
(Ohs Shei dl = 
On 5 - 
De siyA: oe 
O26 a 


DS - Devon Silt - normally consolidated 
DS*- Devon Silt - overconsolidated 


MDS - Modified Devon Silt - overconsolidated 
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Theoretical 

a, ([om/Ymin] 
oul 03 oh D 
07220 226 
0203 208 
0.261 268 
O23 6 242 
0233 239 
0233 239 
Omieo2 239 
G20 232 
0.246 Ze 
0.247 2a 
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each test. An average a value, Qasr could thus be 
calculated. (This method was not applicable to over- 
consolidated samples). A comparison of average a values 


determined in these two ways is presented in Figure 4.6. 
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Comparison of 
Average Frost ; 
0.6 Penetration Rates 
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FIGURE 4.6 


(2) Sample Water Balance 


At any given time during the freezing test 


volume change of the sample is a result of the volume 
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s)he 
change accompanying freezing of the porewater and of 
the volume of water sucked into or expelled from the soil 
sample. It was assumed that the volume change due to the 
freezing of the soil solids is negligible. The volume 
change accompanying freezing of the porewater was assumed 
to be equivalent to the volume of water sucked into or 
expelled from the soil sample, and to the volume displaced 
by the heaving soil surface. By knowing the depth of frost 
penetration and the porosity of the soil, the volume change 
of the porewater (in both frozen and unfrozen states) may 
be equated to (or balanced with) the actual change in sample 
volume. 


As illustrated in Figure 4.7, 


Veeets the volume of frozen soil defined by the 


depth, X, of the 0°C isotherm. 


is the volume portion of Ve that is actually 


frozen at. 0°C. 


V.. is the volume portion of Vie that remains unfrozen 
at nO CC: 

V is defined as a volume of soil that is completely 
PLOzenvate0.ce 

AV. is a volume of water that is sucked into the 

soil sample during the freezing process. 


V is the volume displaced by the heave displacement 


h 


of the soil surface. 
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FIGURE 4.7 


The volume of water in the unfrozen soil zone is assumed 
constant throughout the freezing test and therefore is not 


included in the calculations. Hence we may write: 


m2 
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where Nee assuming Was 0, may also be written 


ee =) 050.9 sn xA ai 


The total volume increase of the soil sample is directly 
measured in terms of vs and AV, - If this volume increase 
is assumed to be equivalent to a known (defined) volume of 


frozen soil we may write: 


O09 ene Vs Vs 2 eer An 4.16 
n ak 


By definition V is equivalent to Ve- Hence by substituting 
V into equation 4.14 the volume of unfrozen water in the 
zone of frozen soil, Mae may be determined. ue is usually 
expressed as a percentage of the weight of soil solids and 
is defined as the unfrozen water content of a soil (Ws). 
Unfrozen water content calculations were made for 
several tests and the results summarized in Table 4.2. 


A sample calculation of Wa is given in Appendix Bb, (Section 2)’. 


(3) Porewater Pressure, Po 


Porewater pressures that theoretically developed 
during various freezing tests were calculated using equation 
2.18.) The calculated Ps were then substituted into 
equation 2.5 and a theoretical K calculated. Results of 
these calculations together with experimental observations 


made during the tests are summarized in Table 4.3. A 
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TABLE 4.2 


Unfrozen Water Content Calculated 


From Sample Water Balance 


Soil Test Unfrozen Water Content 
Thermocouple Measured 
DS E-1-3 hg eke Spyh 
DS E-1-4 85 20 Sei) 
DS* F-1-3 520 3205 
DS* F-2-3 10.00 6.40 
Ds* F-5-3 6.00 ~ 
MDS G-7 14.70 = 
MDS G-8 13.0 = 
DS - Devon Silt - normally consolidated 


DS* - Devon Silt - overconsolidated 
MDS - Modified Devon Silt - Overconsolidated 
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sample calculation of Po and K is presented in Appendix 


Pree CoeCeLOlms) 


(4) Heave Calculation 

Heave was calculated for several freezing tests 
using equation 2.25. Results of these calculations along 
with observed heave values are summarized in Table 4.4. 
A sample calculation of heave is presented in Appendix B, 


(Section 4). 


(5) Heave Pressures 
Results of the two heave pressure tests, Tests 
F-5-5 and G-8, are summarized in Table 4.5. Also listed in 


the table is the volume change due to freezing, V equal 


el 
to the volume of porewater expelled plus the heave 


displacement volume, and the theoretical 9% volume change of 


the irozen sorl zone, V A comparison of V, and Vy indicates 


om F 
the unfrozen water content of the soil, Wie Theoretically 


for pte 0 percent Ve should equal Vo: 
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TABLE 4.4 


Observed and Calculated Heave 


rest Ps 

[Kg/cm?] Observed 
D-1 1.04 0.040 
D-2 1.04 0.040 
E-1-2 0.76 0.452 
E-1-3 P25 O2LZ6 
E-1-4 I= 30 Oe 12 
F-1-3 07.50 0.114 
F-1-4 0.74 0.088 
F-1-5 Ons 0.061 
G-7 2.95 02058 


Ottawa Sand 


Heave [cm] 


Devon Silt - normally consolidated 


Devon Silt - Overconsolidated 


Modified Devon Silt - Overconsolidated 
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CHAPTER 5 
DISCUSSION AND CONCLUSIONS 


5. awascussion 


5.1.A. Equipment 


Once operational, the equipment performed 
Satisfactorily. Reproducibility of test results, as for 
Tests E-1l-1l and E-1-2 and Tests E-2-1 and E-2-2, was good. 

The major aspect of the equipment requiring 
improvement is the pressure sealing system around the 
freezing piston. Some problems occur as a result of the 
sliding O-ring design. One problem is that initial 
pressure cee to the soil sample cannot be maintained 
in the sample container. When pressure is first applied 
to the soil sample a B value of 0.95 is usual but with 
time this pressure dissipates (B = 0.85) and remains at 
this lower pressure. Secondly, whenever pressure is 
applied to the soil sample a small portion of the soil 
sample is extruded upwards past the sliding O-ring. 
Because of this capability of the porewater to extrude 
past the O-ring, some drainage of the top of the sample 
occurs. This water at a later time may go back into the 
soil resulting in a higher water content in the upper areas 


of the soil as compared to the remainder of the sample. 
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This excess water affects the freezing rate of the 
soil and is responsible for the bi-linear penetration of 
the freezing front with respect to the square root of time. 
see Rigurendes theluitiatly the,nate of penetration is 
rather slow, as evidenced by the flat slope of the X, Yt 
curve, whereas below this 'wet' zone the penetration rate 
is nearly twice as fast. The average of the two 
experimental penetration rates, however, compares favorably 
with the theoretical rate of penetration. The difference 
in freezing rates between the upper portion and the 
remainder of the sample does not appear to be reflected in 
the heave or net change in volume porewater results. 

Friction occurred in three areas of the equipment 
but was minimal. The first location was between the 
freezing-piston rod and the O-ring seal. Lubrication was 
frequently applied to this area and any friction that 
developed would be minimal. Friction also developed around 
the sliding O-ring seal located in the freezing piston base 
plate. This area was also lubricated to minimize friction. 
Some friction may also have developed in the sample 
container due to the movement of the soil mass. To minimize 
this friction however, a latex membrane, encasing the soil 
sample, was separated from the sample container by a layer 
of grease. Pressures as low as 2 pounds per square inch 
were 3ufficient to extrude a frozen soil sample (plus 


membran the sample container. 


o0L | . " 
oft Yo siex PRiNse2? od edoStis werew enenie SAAR” 
io aolstsaseneq xeonhi~id- odd xo? eidteaoqaos ai bas Lioe 
-9tt to Jooxr stsupe aid oF toegaex dtiw snoxt patsees? ont 
at noi¢srteneq 30 ster odd YEisisinI .2@.8 sawpkt se2 
+ \x sdt 20 eqola tsii ele} yd bsonebive es ,wole xedts1 | 
ets toitsuteneq sit snos ‘tow' eid? woled eseiedw ,eviIwo_ 
ow odd to apsxove SMT .deed as solwt ylusem al 
yidsrovst astsqmos ,rsvewod ,estsi dott s1tensg fejnemizegxse 
ecmersetiib sAiT .noitsitemeq to etex Iaotisexosds odt Atiw 
eit bas noisiog tsqqL ont nsewtsd estis1 piizesxt at 
ai Setoeltes sd oF tssqqs jon aseb elqmse alt to tobaismss 
-etivess testswsi0q omlov ai spasdo ton to evsesd is 
jnomqiups edt to ese1s sexdt ni baetavo0c noltoixz4 . 
add meewisd asw ncitsool texk? ofl? .Leminim esw sud 
esw aoitsofidul .lees2 pati-O0 sft Sus bor aateiq~patseert 
tend cnortoixi yas Bers s9%6 eidd oF beliqas yltneupett 
bavous beqoleveb oels aoitoixt .Isminim sd bluow beqolevel | 
seed mosata paissext offs ni betsool [ase paisz-O patbile ods 
-MOLI9L7% osimiaim of Boeend seta oeis sew, sexs eldT -stsig 
siqmse edt mi Beqoleveb 'sved oels ysm noitoixt ance 
osimtnim oT .sasm Lioz sit to taomevom sit o¢ sub tenistnoo 
fioe oft eateasone onsidmem xose1 & ,xevewor sofjot1t aids 
teyel 6 Ya tonisinoo slqmse oft moxt bessieges ssw sigqmse | 
font ezsupe 169 ahnbed S es wol as eeuvees1d .sesexp to. 
eviq) slqmse cam Teso1t s Sbustxe ot sar tokue etow 
weniadads sigmse oft nexdmem 


101 
The author believes that errors in experimental 
results due to the pressure seal at the freezing piston 
base plate and due to friction were small and likely did 
not influence the freezing behavior of the sample signifi- 


cantly 


5.1.B. Testing Procedure 

The application and maintainance of a step 
temperature on the soil surface was effectively achieved 
by the constant temperature bath/circulator apparatus used 
in the testing program. This method afforded excellent 
control over the freezing process, though not the freezing 
rate, as recommended by Penner (1972). The freezing rate 
varied from sample to sample depending on the water content, 
density and of course step temperature. The actual step 
temperatures used were chosen arbitrarily on the basis of 
practicality. The experimental step temperatures agreed 
fairly well with temperatures that were reported in the 
literature. Typical temperatures reported in the literature 
were —2 to 12°C by Kaplar (1971) and -17°C by Sutherland 
andsGasking G29 73)-. 

Although the length of the experimental freezing 
tests were short, the tests proved capable of demonstrating 
the freezing behavior of a soil. It was found that as the 
tests progressed the freezing behavior trend that was 


established in the first hour of testing continued for the 
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duration of the test and hence indefinitely depending on 
the sample height. The short term tests were thus 
justified. 

The freezing tests, varied from 120 minutes to 
360 minutes, were usually 160 to 200 minutes in length. 
The length of freezing tests were thus in the order of 1 
hour to 10 or more hours shorter than tests reported in 
the literature by people and agencies such as Penner 
(1972), Williams (1967), and the Corps of Engineers, U.S. 
Army (see Townsend and Csathy 1963). For each freezing 
test the results recorded at 160 minutes (after application 
of the step temperature) were recorded and used for the 
comparison of the various freezing tests. The time, 160 
minutes, was chosen arbitrarily by the author since the 
freezing front was located in the middle third of the 
sample (height) and any effects due to the proximity of 
the bottom or top of the sample were minimal. 

If the tests were not terminated at 160 minutes 
they were terminated when the freezing front was 1.6 
centimeters from the base of the sample. The freezing 
front was not allowed to penetrate any deeper to prevent 
damaging the base transducer. If this depth of penetration 
was reached in less than 160 minutes, as in the case of the 
large step’ temperature, the test results were extrapolated 


toel6éo minutes. 
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As shown on the data summary tables the initial 
water content and void ratio varied for each test. This 
variation occurred as a result of the sample preparation 
procedure. A slurried sample was prepared at a water content 
of 45% (for Devon silt) and then placed in the sample 
container which was partly filled with water. The variations 
in initial water contents listed arises from the different 
amounts of water absorbed by the slurry prior to testing. 
The author feels that these variations in water content 
and void ratio did not affect the freezing test results 
Significantly and is only of interest in terms of sample 


preparatiomicontrol: 


yee. Experimental Results 


Experimental results indicate that the freezing 
behavior of a soil can be described in terms of net change 
in sample porewater volume, porewater pressure and 
heave. Results also indicate that soil type, grain size, 
grain size distribution, permeability, stress history, 
overburden pressure and step temperature influence the 
freezing behavior of a soil. 

The difference in the freezing behavior of coarse 
and fine-grained soils is clearly shown in Figure 4.1. 

The Ottawa sand data indicate that under all freezing test 
conditions porewater is expelled. The freezing behavior of 


Devon silt however depended on the test conditions. The 
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Devon silt expelled water under certain effective stress 
conditions, roughly 15 to 20% (by volume) of that of the 
Ottawa sand. This difference in freezing behavior can 
be attributed to differences in shape and surface area of 
the soil particles, and permeability as determined by 
grain size and grain size distribution. 

The freezing behavior of Devon silt shown in Figure 
4.1 should be noted. The results show that the freezing 
behavior of Devon silt is determined primarily by the grain 
size, as implied by the capillary theory, but is strongly 
dependent on the applied effective stress during freezing. 
As illustrated in Figure 4.1 Devon silt may behave in three 
different ways during freezing tests: (1) increase total 
sample porewater volume, (2) decrease total sample pore- 
water volume, or (3) maintain the total sample porewater 
volume. These freeze test data further imply that a 
linear relation between change in total sample porewater 
volume and effective stress exists, and that the freezing 
behavior is a continuous phenomenon with a smooth transition 
from one behavior mode to another. 

The results show that at a particular pressure, Por 
there is no net change in total volume of sample porewater 
Under these freezing conditions it is likely that all volume 
change, due to the 9 per cent volume increase at phase 


transformation, be expressed in terms of heave by the soil 
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sample. At pressures less than Por water is sucked into 
the sample by the freezing front, increasing the total 
sample porewater volume. This action, if allowed to 
continue under the same conditions, could lead to ice 
segregation and lensing. At pressures greater than Por 
porewater is expelled from the sample thereby decreasing 
the total porewater volume. Under these conditions in 
situ freezing results; ice segregation and lensing can 
eccury 

The preceding results are consistent with the 
Theory of Freezing Soils proposed in Section 2.2 and 
predicted by equation 2.5. When no net change in sample 
porewater! occurs, the porewater pressure Pou is equal 
€O9zero, Or Ps of equation 2.5 and of the freezing sample 
is equal to zero. The effective pressure, Por of Figure 
4.1 then is equal to P, which in this case is also equal 
EQOCKEOLrecquationn2 .SoaaThusetermseot Keand Py are inter- 


changeable. 
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At Po = 0 equation 2.5 and the accompanying theory 


predict that water would be sucked into the soil sample 
and ice lensing would result. This phenomenon was not 
verified experimentally and warrants further study. 

It should be noted that in all tests shown in 
Figure 4.1 and in the testing program the freezing plane 


was continually moving through the soil. At first this 


te. ae 

; ie 
odnt bexoue et 9tew 4% asdd asel esmvasorg 26. -eLqmas 
Letos ontt priasoroni anaes paiseost eds yd eigmss. oft 
ot fewolls 24 ,aottos eidT? .emulov 1stewst0q elqmss 
aot ot Deel Bbiyop ,a@acisibaos smeae oft isbee sumitmoo 

1,7 “edt aetsonp asauapeng tA -paiedel brs noitspsipes 
pitasexceb ydexedt+ olqmke ont mot belleqxe at 167 BWS10g 
nt anoitibaod aeons xoball somiioy iedsweneg {stot ens 
nego patesel bas noitspewese soi ;ativeset pnissest utte 

| | | «1U090 
efit ritiw tnevatenoo ots edivesx paibsoeag ect 
bee SS noitoee ak beeogotq elio® paises to yxzoedT 
elqmse ai opasdo ton on medW .2.8 moitsups yd hetpibeag 
lavps, at vag Ssive2esiq itstswetog eft ,erveDpe iIsedswet0g 
elqmsa putssext odd 26 bas 2,.£ sottsups to wt 7 ,Or9s. 02 


sxvpid toe % ,etweesiq svisostic siT .O1S8 od Lseyps. el. 


‘oO 
leupe oAls ai senso airs m4 toi clw 33 os iapee ei nena Ld 
-asdénl 926 oo Sas A to eurtet auiT .2.8 moitaups to A os 
| | midnemenie | 

ytood? priynsqmonss eid bas @.% aoivaups.0 = 4 JA 
slamee {toe edt otnl bedows asd bivew t]etaw gadt solbexrg 
al eaw nonemomariq eitiT .tivees —— eaenaes ssi Bas 
ybose sortaud. atasi1ew bas es kasnomhxeaxe beitiisv 

mui awode ateot Ife wt ted? beton ed blyona tI 
esiq paisesit odd mexpoxq puitees odd ai bas 1. oxupht 
eins taxi? aA -tioe eit dgworit paivom uitsunidaem caw 


ae ei bet ect ‘1% 


106 
appears inconsistent with the theory and observations 
reported in the literature that ice lensing or the sucking 
in of water during freezing occurs only when the freezing 
front is stationary. However, the apparent inconsistency 
is likely due only to the scale of the sample height, step 
temperature and freezing time used in the freezing tests. 
The freezing front in these tests is advancing because the 
rate of heat removal is greater than the rate of heat 
supply. For a similar sample of infinite height and 
subjected to the same step temperature, the freezing front 
would advance and suck in water until the heat exchange and/or 
pressure conditions changed significantly. This freezing 
phenomenon is de ained as ‘imperfect segregation' by Arakawa 
(1966). As the freezing front advanced the rate of heat 
removal would decrease in proportion to the depth until the 
rate of heat removal was equal to the rate of heat supplied. 
The freezing front would then become stationary. At any 
time during this advance of the freezing front the sample 
would continue to suck in water as long as the overburden 
pressure P. was less than Po At this time the freezing 
front would become stationary and continue to suck in water, 
assuming the difference in pressure P; = Pe of equation 2.5 
was still less than (or equal) to K. 

Equating equations 2.2 and 2.5 we see that K may be 


calculated in terms of: 
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Sutherland and Gaskin (1973) suggest using: 


a) 


ae = 0.03 N/M = 3.06 x 10 Kg/cem 

r/x; = 5.4 562 
where ee ~ stress across the ice-water interface 

a ~ radius of the ice crystal interface 

1G - radius of the soil pore = 0.5 (Dj 9) 


From Appendix A, D of Devon silt is 0.00016 mn. 


10 
Subst LUUCVONMOl rei nUuO venuatvOue o.2,) aC eduatloN sage. nto 
equation 5.1a value of 41 kg /om* is determined for k. 

The experimental K value for normally consolidated Devon 
Ste er bre o0 kg/em*. 

Further to Figure 4.1 it should be noted that the 
linear relation between net change in sample porewater 
volume, AV, and effective stress, o', exhibited by Devon 
silt also holds true for the finer-grained modified Devon 
Silt. The slope of the AV - oc! curve of the modified Devon 
silt is less than, and the effective stress intercept, Por 
is greater than the corresponding values for Devon silt. 

It appears then from Figure 4.1 that the slope of the AV - o' 
curve and the es intercept are unique characteristics of a 
soil type. It would also follow that the characteristics of 
the AV - o' plot are dependent on soil type, grain size, 


grain size distribution and percentage of fines. Again 
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the role of soil permeability should be noted. 

It is also evident from the test data that the 
freezing behavior of a soil is dependent on the stress 
history of that soil. From Figure 4.3 it can be seen that 
the slope of the AV - o' curve increases, and the aes intercept 
decreases as a result of overconsolidation of the Devon 
silt. The plot also shows that the slope of the AV - o' 
curve increases and the ae intercept decreases with an 
increasing degree of consolidation. The AV - o' curve 
plotted for Test Series F-5 where the preconsolidation 
pressure, Bey equals 8.79 Kg/em* is inconsistent with the 
above observations. This inconsistency is likely due to 
experimental error. 

Replotting the results shown on Figure 4.3 indicates 
that ee is inversely proportional to the log, 9 of the ratio 
of B/ee where Po is the preconsolidation pressure. It 
can be seen from Figure 5.1 that once the a= of a soil 
type is known the influence of stress history on that Po 
can be determined. Results gained from tests run on Devon 
silt at an increased step temperature and tests run on 
modified Devon silt do not correspond to the Devon silt 
data points. This would indicate that the fa log, 9 Po’ s 
curve is unique for a soil type being frozen under certain 
temperature freezing conditions. Figure 5.1 has many 
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Further to Figure 4.3 the effects of the 
difference in testing procedure should be examined. For 
the test series run on normally consolidated samples a 
new sample was used for each freezing test. For the series 
run on the overconsolidated samples, one sample was used for 
each series, and therefore subjected to several freeze-thaw 
cycles. As shown by the test result curves the AV - o' 
relationship is linear once the sample had been subjected 
to one freeze-thaw cycle. However, the slope of the same 
curve is significantly steeper between the lst and 2nd 
freeze-thaw cycles. This change in slope indicates that the 
effect of the freeze-thaw consolidation of a soil is much 
greater for the first freeze-thaw cycle than for subsequent 
cycles. 

Hence a comparison of the es of normally consolidated 
samples determined on the initial freezing of a sample 
should not be stricly made with the Po of an overconsolidated 
sample subjected to several freeze-thaw cycles. Consequently 
the normally consolidated samples should have been subjected 
to at least one freeze-thaw cycle prior to determining a ays 
value for comparison with the Po of overconsolidated samples. 
The author recognizes this shortcoming in the results. The 
author also believes that if the more valid procedure for 
normally consolidated soils had been followed the result 


would have been a Po intercept occurring at a larger 
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effective stress and a flatter slope on the AV - o' curve. 
This prediction is based on the shape of the lower portion 
as compared to the initial portion of the overconsolidated 
AV - o' curves. The comparison of normally and over- 
consolidated test result data shown in Figure 4.3 is nonethe- 
less useful. 

Freezing tests conducted under close system 
drainage conditions were run to determine if the freezing 
behavior of a soil was affected by the drainage boundary 
conditions. Closed system freezing tests are also a 
practical alternate testing procedure. Test data for 
normally consolidated and overconsolidated Devon silt 
samples are shown in Figure 4.2. Test results for the 
normally consolidated samples are somewhat erratic and 
hence may be questioned. Data for the overconsolidated 
samples are more consistent. The difference between the 
quality of results is believed due to the higher compress- 
ibility of the normally consolidated samples. 

The trends of the closed system freeze results are 
consistent with open system freeze tests discussed earlier. 
The results show that porewater pressures are equal to zero 
at a particular effective stress, Rag at stresses greater 
than Por the porewater pressure is positive (or porewater 
would be expelled from the sample); and at stresses less 


than Py the porewater pressures are negative (or water 
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would be sucked into the soil sample). The Ey values of 
closed system freezing tests do not favorably compare with 
those of open system freezing tests. Table 5.1 summarizes 


the experimental Po values. 


TABLE, 5 2 


Experimental Po Values 


Stress History Po for Freezing Test [Kg/em“] 
Open Drainage Closed Drainage 

Normally 

Consolidated 23 1S 746) 

Over 

Consolidated 

Sea = 7.70 Kg/cm Ont OF 5S 


As can be seen in Table 5.1 the freezing behavior 
of Devon silt appears to be dependent on the drainage 
boundary conditions of the freezing test. 

DUPanGe cher progress: Of lests@i—s-l “and E=3-2mMthe 
drainage boundary condition was changed to further check the 
influence of the drainage boundary condition on the freezing 
behavior. Results of these tests were inconclusive. The 


freezing behavior of the Devon silt was unaffected by the 
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ab 
change in drainage boundary conditions in Test E-3-2. 
However, the freezing behavior of Devon silt in Test E-3-1 
was completely reversed. In Test E-3-l the sample expelled 
water under open system drainage, showed a steadily 
decreasing porewater pressure under closed system drainage, 
and then sucked in water when the drainage system was 
opened again. 

Because of the inconclusiveness of the closed 
drainage system freezing tests, this type of testing was 
abandoned and efforts were concentrated on open drainage 
freeze tests. Further work on closed system freezing tests 
is warranted. 

As winnie os in Figure 4.4 the freezing behavior 
of Devon silt is dependent on the step temperature. When 
the step temperature was doubled the same freezing behavior, 
as discussed earlier,was followed. The slope of AV - o' 
curve was steeper, however, and Be was less than the corre- 
sponding curve for a step temperature of -5°C. 

Devon silt does not completely freeze at 0°C. The 
unfrozen water content of Devon silt was determined in two 
ways: by measurement of the advance of the 0°C isotherm 
using thermocouples and by vernier, and secondly by determining 
a sample porewater balance and back calculating the unfrozen 
water content. Calculated unfrozen water contents ranged 


from 4 to 6 percent. These values compared favorably with 
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an unfrozen water content of 5 percent for a comparable 
silt at the same freezing or step temperature by Anderson, 
Tice and McKim (1973). Because of this unfrozen water 
the total change in sample volume during freezing was 
always less than the 9 percent volume change of water at 
nucleation. 

Based on the general agreement between predicted and 
observed freezing behavior, the porewater pressure at the 
freezing front can be estimated from the net supply and 
removal of heat from thesoil sample. All calculated pore= 
water pressures listed in Table 4.3 were determined by 
equation 2.18. —Substition of the calculated porewater 
pressures into equation 2.5 enabled the author to predict 
freezing behavior. The results of these calculations indicate 
that equations 2.5 and 2.16 are useful. These equations can 
predict the freezing behavior of a soil if the step 
temperature, water content, permeability, thermal properties, 
dry density, water table conditions and the K value of the 
soil are known. 

With the exception of two tests all samples 'heaved' 
during the freezing process. The two samples that did not 
heave were tested under restrained heave conditions. A 
small amount of vertical displacement in these two tests 
did occur however due to compliance. The vertical 


displacement or 'heave' measured during the remainder of the 
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freezing tests was likely due to the 9 percent volume 
increase at phase transformation. Tests in which water 
was being sucked in would develop segregated ice and ice- 
lenses if the tests were allowed to continue for a sufficient 
period of time. The vertical displacement in this case would 
correctly be termed heave as is done cn the literature. The 
samples which expelled porewater during freezing tests 
would displace vertically little more than recorded. These 
samples would not heave. As expected the greatest displace- 
ment occurred under the smallest overburden pressures. These 
results indicate that heave is merely a volumetric response 
along with the expulsion /attraction of water to the freezing 
of the soil porewater in conjunction with the applied 
boundary conditions. The results imply that the magnitude 
of heave is not unique to the particular soil sample 
tested. 

Calculation of vertical displacement, based on 
the volume increase at phase transformation of the amount 
of porewater frozen during the freezing process, was 
attempted. Calculated results were in poor agreement with 
the observed heave. The reason for this discrepency is 
likely because the overburden loading was not. accounted for 
in the calculations. 

The 'heave' pressure results are summarized in 


Table 4.5. The summary shows that both samples expelled 
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water during the freezing tests or in other words neither 
sample was developing ice lenses during the freezing tests. 
Therefore the measured pressures were a result of the 9 
percent volume increase at phase transformation. For the 
sake of clarity these pressures will be referred to as 
transformation pressures. 

Since the applied overburden pressures used in 
both tests were such that no expulsion of porewater would 
occur during freezing, the transformation pressures were 
responsible for the expulsion of porewater This can be 
seen from the test results; the largest volume of porewater 
expelled corresponded with the largest transformation pressure 
measured. The results also show that the largest trans- 
formation pressure was developed by the sample with the 
smallest overburden pressure. Therefore heave pressure 
depends on overburden pressure. Contrary to Gold (1957), the 
‘coarser' Devon silt sample generated larger transformation 
pressures than the finer modified Devon silt sample. This is 
likely due to a higher unfrozen water content and the lower 
permeability of the modified Devon silt. 

Movement of the freezing piston was also recorded 
and summarized in Table 4.5. . This vertical movement is a 
measure of the compliance of the rigid frame. Any upward 
or compliant movement of the rigid frame resulted in lower 


transformation pressures being measured. This reduction 
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in pressure undoubtably happened in both tests but the 
effect of this reduction appears to be minor in comparison 
to any change in overburden pressure. 

In summary, heave pressure is pressure developed when 
vertical displacement is restrained due to the generation of 
transformation pressures plus pressures resulting from the 
development of ice lenses (also called heave pressure). 

From the test results quoted it appears that heave pressure 
is dependent on soil type, permeability, grain size, and grain 
size distribution, and overburden pressure. 

Based on the above summary, results, and discussion 
the following theory on the mechanism of heave pressure is 
outlined. 

Whenever the soil surface is restrained during an 
open system drainage test the ice crystals displace the 
soil porewater. These ice crystals may occur as an ice lens 
or individually as after crystallization. The resistance 
to flow that the porewater encounters is the basis of heave 
pressure. Bounds may then be set on the flow resistance, in 
terms of pressure, of various soils: (1) in coarse grained 
soils with a high permeability the resistance to flow would 
be minimal, (2) in a fine grained soil with a very low 
permeability near infinite pressures are required to initiate 
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Therefore in the case of the granular soil, as the 
soil freezes the ice would fill the voids and the displaced 
porewater would be able to move freely throughout and even 
leave the soil sample. The pressures generated on the soil 
surface would be negligible as would be the vertical 
displacement. For the fine grained soil however, the pore- 
water being displaced by the growing ice crystals (or lenses) 
would encounter significant resistance to flow. Initially 
no flow would occur but as the ice crystals enlarged increased 
pressures in the voids would result. When the pressure in the 
void became larger than the soils' resistance to flow the 
porewater would move out of that pore space. The pressure 
generated to cause this flow could be measured on the soil 
surface and would be defined as heave pressure. 

In summary, heave pressure would be measure of the 
degree of the free drainage capability of a soil system. 

In the field then if a soil permitted free drainage, open 
system drainage conditions would be approached; whereas if 
the soil permitted practically no flow then closed system 


drainage conditions would be approached. 


Soe meCONGLUSIONS 

Experimental results indicate that the short term 
freezing tests conducted during the testing program proved 
to be an effective procedure in determining freezing behavior 


of a soil. Therefore the running of prolonged freezing does 
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not appear necessary. The application and maintainance of 
a step temperature on the soil surface, as opposed to varying 
the surface temperature to maintain a constant rate of 
penetration of the freezing front, is a very effective method 
of maintaining experimental control over freezing tests. 
Application of a step temperature is also a simpler and more 
practical method of freezing samples. 

For the soils tested, the running of open system 
drainage tests, measuring net change in sample porewater 
volume and heave, and closed system drainage tests measuring 
porewater pressure and heave, proved to be effective methods 
of determining the freezing behavior of a soil. The 
diversity of experimental results obtained and of testing 
boundary conditions used show that the freezing behavior of 
a soil cannot be accurately predicted on the basis of one 
soil property, i.e. grain size, or one freezing characteristic, 
i.e. heave pressure. Therefore in describing the freezing 
behavior of a soil the boundary conditions such as step 
temperature, drainage and pressure conditions, and soil 
properties should be listed. Care should be taken when 
comparing data obtained from open and closed system freezing 
tests as experimental results indicate the freezing behavior 
of a soil is influenced by the testing boundary conditions. 
Further study of these affects is needed. 

Soil properties that affect the freezing behavior of 


a soil are: soil type (specific surface area), grain size, 
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grain size distribution, density, porosity, water content 
and permeability. Test results also demonstrated that the 
freezing behavior of a soil is dependent on external factors 
such as: overburden pressure, stress history, drainage 
boundary conditions and step temperature. The relative 
importance of these factors for the soil types tested is as 
follows: 

Cohesionless soils used with no fines, expel water 
upon freezing. This phenomenon appears insensitive to 
changes in overburden pressure, porosity or dry density. 

The freezing behavior of a fine grained soil (silt) 
is not only dependent on soil properties but is also 
dependent on external factors. The freezing behavior 
exhibited by the silt varied from sucking in to expulsion 
of porewater in open system freezing tests. 

A definite relationship was shown to exist, for the 
soils tested, between the net change in sample porewater 
volume, AV, and effective stress, o'. At a particular stress 
level, PS! no net change in the sample porewater volume 
occurred. (This however was not verified experimentally). At 
stresses less than Po water is sucked into the sample which 
would result in ice segregation, lensing and significant heaving 
if allowed to continue indefinitely. At stresses greater than 
P_, porewater is expelled during freezing resulting in situ 
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freezing and insignificant heaving if allowed to continue 
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indefinitely. The test results indicated that the freezing 
behavior of the Devon silt was’a continuous phenomenon with a 
smooth transition from one type of behavior to another. 

Results of tests run on the modified silt further 
indicated that the AV - o'relation was unique to a 
particular soil type and that the slope of the curve and the 
Ee intercept is influenced by grain size distribution and 
permeability. It was shown that the slope of the AV - a' 
curve decreases and that the By intercept increases as the 
percentage of fines in a soil increases, or in other words 
as the permeability of a soil decreases. 

Besides being stress dependent the freezing behavior 
of the silt was also shown to be dependent on stress history. 
Overconsolidation of the silt increased the slope of the 
AV -="o.""curve ‘and. decreased ‘the Bt intercept... Furthermore 
it was shown that the slope of the AV - o' increased and the 
Py intercept decreased as the degree of overconsolidation 
increased. 

When the 'stress history" data were replotted it was 
found that a was inversely proportional to the logs (P /Po) 
where Po is the preconsolidation pressure. The data indicated 
that this curve was also unique for a soil type and was 
dependent on step temperature. This relation between Ps and 
the log, 9 (P /P,) has important practical implications, ive. 
in estimating the depth of burial (overburden) required to 


prevent frost heaving of a pipeline. 
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Test result data also demonstrate that the freezing 
behavior of a fine grained soil is dependent on step 
temperature. It is shown that the slope of the AV - o' 
curve increased and the Pe intercept decreased when the step 
temperature was increased. 

Closed system freezing tests on the silt demonstrated 
that freezing behavior could also be monitored in terms of 
porewater pressures at the freezing front. Freezing behavior 
of the silt exhibited the same trends but did not duplicate 
the behavior shown in the open system tests. Based on the 
experimental data it appears that freezing conditions are 
not identical under open and closed drainage conditions. 

The validity a6 the assumptions and theory of 
freezing soils described earlier was demonstrated by the 
experimental results. Theoretical predictions of freezing 
soil behavior developed in terms of the relative heat flow 
into and out of the soil sample and of the stress difference 
across the ice-water interface were verified by experimental 
results. 

Little success was achieved in predicting the heave 
of a soil. This is likely due to the limitations of the 
theory and due to the influence of the unfrozen water content 
and the omission of overburden stress in the development of 
the heave equations. The test results indicated that heave 
was not a property of a freezing soil but is a response of 


a freezing soil dependent on the overburden pressure and 
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drainage boundary conditions. 

Similarly the test data showed that heave pressure 
is not a soil property but is also a response to the 
restrained displacement, drainage and overburden boundary 
conditions. The amount of porewater expelled during the 
heave pressure tests appeared to be linked to the heave 
pressure measured. 

The experimental data demonstrated the diversity of 
freezing behaviors that a soil (fine grained) can exhibit 
depending on the boundary conditions. The experimental data 
also showed that for the soils used the freezing behavior of 
granular soils does not vary within the same limits as the 
fine grained soils. This evidence clearly indicates that the 
present frost susceptibility criteria, based on grain size, 
is inadequate and needs revision. It was clearly demonstrated 
that a fine grained soil cannot be (in absolute terms) classed 
as frost or non-frost susceptible. The freezing behavior of 
a soil is dependent on many factors and in order to classify 
the degree of frost susceptibility of a soil these factors 
must be dealt with. A new frost susceptibility criteria 
should set definite limits on freezing soil behavior in terms 
of soil properties, drainage, temperature and applied stress 


boundary conditions. 
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Le ample Calculation -—- Theoretical Frost Penetration 
Rate 
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By using Simpson's Integration formula the 
T~time plot shown in Appendix E, (Figure E-1-2) was 
integrated and an equivalent temperature T, was determined 
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2. Water Balance = Sample Calculation 


Consider Test E-1-4 
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By using Simpson's Integration formula the Ts — time 
plot shown in Appendix F, (Figure F-1-4) was integrated 
and an equivalent temperature T. was found or: 
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CONSTANT TEMPERATURE BATH/CIRCULATOR MODEL 334 


Available from HOTPACK (CANADA) LTD. 


385 Phillip Street Nw, Waterloo, Ontario. 


TEMPERATURE RECORDER 


HONEYWELL ELECTRONIK 15 STRIP CHART MULTIPOINT 
RECORDER 15303836 -24-04-0-000-004 10 
Available from HONEYWELL INDUSTRIAL PRODUCTS GROUP 


Wayne and Windrim Avenue, Philadelphia 44, Pa. 


"BELLOFRAM' FLEX-WEAVE MEMBRANE 


Catalogue No. 3-400-244 DBT 
Available from ANTEUS LABORATORY EQUIPMENT INC. 


Malden Mass. 02148 


JAYCOR UNSATURATED POLYVINYL CHLORIDE (UPVC) PIPE 


Thermal conductivity at 20°C = 4 x Nop: [ cot 


Specific Heat 0.24 Caan 


Specific Gravity 1.47 [gm/em?] 


Modulus of Elasticity - tension 350,000 [p.s.i.] 


- flexure 400,000 [p.s.i.] 


144 


.olrsinO (,oolxeteW ,.u teers2 gillidd ese 


THIOTLTIUM TAAHD PIATS 2@L ALMOATOAIA LIAWYSMOR 
Of BO0~000-0-80-BS- DEGEOEZL AAGHOOSA 


- 
a 


qUOSD 2TOUdGOST JALATAUGMI GISWYSMOH mos? sidslisvA 
-59 ,68 sidqlobslidd ,sunevA miabnkw bas omyaw 


a 
et Fee 


= 7 

a 

SMARGMAM TVASW~XaIG 'MANIOIIEa* : 

- 

Pad BbS-008-E .o% supolsisd 7 

.OMl TMAMTIVUON YSOTAHOGAI 2USTUA mox? oldslisvA | 
SbISO .2esM neblsm : 


= = 
@IT¢ (OVI) SaLAOUHD IYUIVYAOd aaTARUTAeMU RODYAL 


iss, %- > 
as Of xb’ = D°OS 45 viivlseabaco Lamsedt ‘7. 


igotge] 88.0 s50H ottivege 
( mo\me! Ta.L ydivserd oliioege 


[-i.8.q) 000,G2€ motenss - ytiotzesia lo aulubom 
{.L.a.q] 000,008 exuxel> = 


145 


‘ -5 ,in 
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41°F to 140°F 


Available from JOHNSTON INDUSTRIAL PLASTICS 
9537 - 62nd Ave. 
Edmonton, Alberta. 
NULLMATIC PRESSURE REGULATORS MODEL 40 


Available from MOORE PRODUCTS,. Rexdale, Ontario. 
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